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SUMMARY 


Experiments of transonic flow past a circular arc profile show 
that the shoek-wave pattern and the pressure distribution are 
strongly dependent upon the state of the boundary layer. A 
change from laminar to turbulent boundary layer at a given Mach 
Number changes the flow pattern considerably. 

Shock waves can interact with the boundary layer in a manner 
similar to a reflection from a free jet boundary. These shock 
waves are not distinctly discernible from pressure distribution 


measurements. 


INTRODUCTION 


MM EXPERIMENTAL INVESTIGATIONS of transonic 
flow phenomena have been based on the assump- 
tion that the Mach Number is the only important param- 
eter. Reynolds Number effects were assumed, in 
general, to be negligible. There are,. however, certain 
reasons to suspect that Reynolds Number effects can 
be appreciable and that the boundary layer does in- 
fluence the transonic flow pattern. One of these rea- 
sons is that potential theory has shown that smooth 
transonic (that is, subsonic-supersonic-subsonic) flow 
is possible with local Mach Numbers well in excess of 
unity (for example, reference 1). Experiments? to date 
have failed to show the existence of large local super- 
sonic regions without the occurrence of shock waves. 
This discrepancy betweén theory and experiment can 
be due to the fact that the potential solution for shock- 
free transonic flow is a possible but unstable solution, 
while the flow with shock wave is more stable and occurs 
in experiment; or the discrepancy may be due to the 
fact that certain basic assumptions of the theory (i.e., 
steady, nonviscous flow of a perfect gas) are not met in 


the experiments. Of these assumptions, the one that 


| is most likely to fail is that of zero viscosity. In nearly 
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all cases, transonic flow occurs near a solid boundary. 
Thus, the supersonic zone will be in contact with a 
laminar or turbulent boundary layer, the existence of 
which is neglected in potential theory. 

These two possibilities may be illustrated with well- 
known examples from real fluids. In pipe flow laminar 
motion is a possible solution for all Reynolds Numbers. 
Actually, beyond a certain Reynolds Number turbu- 
lent flow is always found experimentally. This is 
due to the fact that laminar flow is not a stable solu- 
tion. 

A different example is flow in a divergent channel or 
diffuser. It is well known that the flow pattern for a 
nonviscous fluid and for a viscous fluid are radically dif- 
ferent, even for vanishingly small viscosity. The actual 
flow pattern of the real fluid cannot be approximated 
by perfect fluid flow. 


In fact, it is fortunate that in most subsonic aero- 
dynamic problems full use can be made of the classical 
boundary layer theory, and that pressure forces can be 
computed on the basis of the perfect fluid flow, i.e., 
potential theory. The influence of viscosity and tur- 
bulent shear is generally restricted to a narrow zone, the 
boundary layer, and is important for friction effects 
only. It will be seen that in some problems of transonic 
flow, and probably for some supersonic flow problems, 
the classical boundary layer assumptions fail so that 
results derived from boundary layer theory become 
questionable. 

The investigations reported here have as their ulti- 
mate goal the understanding of the formation of shock 
waves in transonic flow. The reasoning outlined above 
led to the belief that boundary layer effects are of para- 
flow phenomena. 

boundary 


mount importance in transonic 
Hence, the interaction between 
and shock wave was the object of the first investiga- 


tions. 


layer 
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Fic. 1. The GALCIT 2- by 20-in. transonic wind tunnel. 


The work presented here was carried out during 1945* 
at the California Institute of Technology under the 
sponsorship of: the Air Technical Service Command, 
Wright Field. 


APPARATUS AND METHODS 


Wind Tunnel 

The 2- by 20-in. transonic wind tunnel, built speci- 
fically for investigations of this type, was used in these 
experiments. A general layout of the tunnel is shown 
in Fig. 1. The tunnel is driven by two turboblowers 
that can be operated either in series or in parallel. The 
power is sufficient to choke the tunnel with any model. 

The emphasis in the design of the tunnel was on flexi- 
bility and quality of the airflow, i.e., low turbulence 
and steady and easily adjustable flow. Considerations 
of energy ratio were of secondary importance. Hence, 
screens are used freely in the circuit to reduce the turbu- 
lence level. The speed is regulated by throttling and 
by-passing, the blowers being driven by constant speed 
induction motors. An extremely fine speed control 16.2. The test section of the GALCIT 2-by 20-in. transonic 

. . i . ‘ d wind tunnel with the 3-in. by 12 per cent biconvex airfoil in place. 

between VM, = 1 and WM, = 0.1 is obtained in this man- “ 





* The research was actually initiated in 1944. The first year 


was spent in the design and construction of the wind tunnel and €T- Details of the tunnel construction can be found in 


reference 3. 


Fic. 4 


equipment. 
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Vertical pressure distribution on tunnel wall empty 
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Test section wall showing pin vises (1) and phosphor bronze wires (2) used in mounting the model (3). 


Test Section and Model 


The test section is seen in Fig. 2. The dimensions 
are 2 by 20 by 36in. The section widens from 2.00 in. 
at the entrance to 2.21 in. at the exit to keep the Mach 
Number distribution uniform. The maximum varia- 
tion in Mach Number along the section at M = 0.85 
was found to be 2.5 per cent. 

The model support is seen in Fig. 3. The airfoil is 
held between the glass walls by friction. Two fine 
wires stretch upstream along the walls of the test sec- 
tion. These wires are for safety only and do not carry 
much load. The wires are in the wall boundary layer 
and do not, therefore, disturb the flow past the model. 
This type of model support was found to be suitable for 
the present experiments. The drilling of holes in the 
glass walls is avoided, and the position of the model in 
the test section can be changed easily. 


Humidity Control 

For the investigations presented here it was found 
that raising the stagnation temperature in the tunnel 
gave adequate humidity control. The cooling water 
for the radiator can be regulated easily and the temper- 
ature in the tunnel adjusted to within 0.2°. The 
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measurements were all made at a stagnation tempera- 


ture. 
to = 55.5°C. = 0.2° 
The dew-point temperature was always found to be 
Jv, = 14°C. = 2° 


The stagnation temperature was measured by means of 
a resistance thermometer in the settling chamber. The 
saturation temperature was measured by means of a 
dew-point meter of the type described in reference 4. 


Determination of Mach Number 


The Mach Number 7 in the test section was ob- 
tained from readings of the difference between the pres- 
sure in the settling chamber fp and the average pressure 
pr in the test section upstream of the model (fo — pr) 
and of the absolute value of pr. My is then given by 
the well-known relation 


po/br = {1 + [y — 1)/2] Mr} O-” 


y = Cp, Cy was always taken as 1.400. 
The pressure distribution normal to the direction of 
flow in the reference plane is shown in Fig. 4. 


Schlieren System 


The Schlieren system employed in these investiga- 
tions is of the two-mirror type. The actual light path 
is complicated because of lack of space. The light path 
The quality of the pictures is not 


is shown in Fig. 5. 
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Schematic drawing of the Schlieren optical system of the 2- by 20-in. transonic wind tunnel. 


adversely affected by the rather complicated setup. The 
photographs are sharp enough to allow enlargement to 
more than six times the original size of 4 by 5 in. 

For continuous observation, a mercury vapor lamp 
(G-E H4) was used. The photographs were made by 
means of a spark discharge of 4 K 10~® sec. duration. 


Pressure Distribution Measurements 


For pressure distribution measurements, one of the 
circular glass plates is replaced by a stainless-steel plate 
upon which the model is mounted. The pressure tubes 
are taken out through the steel side wall. A standard 
multiple mercury manometer is employed for the pres- 
sure readings. 


Laminar and Turbulent Boundary Layer 


Schlieren photographs taken with the knife edge par- 
allel to the mean flow direction—i.e., so that density 
gradients normal to the direction of flow are recorded 
clearly show the boundary layer. In fact, it is found 
that laminar and turbulent boundary layers are easily 
discernible and the point of transition is marked. In 
the laminar case the edge of the boundary layer is well 
defined, whereas in the turbulent case the edge fades 
out into the free stream. This is due to two facts: (1) 
The laminar boundary layer profiles (velocity, density, 
or temperature) have a much larger curvature near the 
edge of the layer than do the turbulent profiles; (2) 
the boundary layer fluctuations in the laminar case are 
of a much larger scale than those in turbulent layers. 
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Hence, the optical “‘surface’’ formed by the boundary 
layer is ‘“‘wavy”’ in the laminar and “‘rough”’ in the tur- 
bulent case. 

The transition point of the airfoil used in the present 
investigation was, except for the very highest Mach 
Numbers, downstream of the shock wave. Transition 
was found to occur at about 80 per cent chord at a Mach 
Number M;, = 0.795 and at a Reynolds Number, 
based on the airfoil chord, of R = 850,000. 

To investigate cases in which the boundary layer is 
turbulent at the position of the shock wave, the follow- 
ing simple procedure was used: A thin ('/;¢ in.), 
pointed flat plate was fastened to the leading edge of 
the airfoil as shown in Fig. 6. This plate, at zero angle 
of attack, does not appreciably affect the pressure dis- 
tribution. It does, however, produce an effectively 
larger Reynolds Number and a turbulent boundary 
layer. Using this scheme, measurements at the same 
Mach Number with both laminar and turbulent bound- 
ary layer could be made and compared. Other obvious 
means to render the boundary turbulent, such as sur- 
face roughness elements and increased tunnel turbu- 
lence, were also used occasionally. 


PURPOSE OF THE EXPERIMENTS 


The experiments were undertaken to study the mech- 
anism of the boundary-layer-shock-wave interaction 
rather than to test specific airfoil shapes. The choice of 
model, therefore, was dictated only by considerations 
of simplicity. A 12 per cent thick, double circular arc 
section of 3-in. chord was chosen, and all experiments 
reported here were carried out with this model at zero 
The critical Mach Number for this 
airfoil is MW. = 0.79. The tunnel was found to choke 
at Mach Number M; = 0.91. The investigations of 
transonic flow thus covered the range between these 
limits. ; 

Wind-tunnel wall corrections are of minor importance 
in the present investigation, since no attempt to simu- 
late free flight conditions is made. The Mach Numbers 
given in the figures are thus the uncorrected Mach Num- 


angle of attack. 


bers. 

Two types of measurements were made, Schlieren 
photographs and pressure distribution surveys. It will 
be seen that, in spite of many prejudices, Schlieren 
photographs can give a clear picture of the flow pattern 
and, furthermore, when properly evaluated, can furnish 
a great deal of quantitative information. Schlieren 
photographs have the great advantage of being the re- 
sult of a nearly instantaneous observation. This is 
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Three-in. by 12 per cent circular arc airfoil + 3- by '/,.-in. flat plate 


important for boundary layer effects, which by nature 
are quasi-stationary and fluctuating. 

Pressure distribution measurements, especially with 
models of small scale, have a large time lag and furnish 
a picture of mean effects only. The difference in time 
lag must be taken into account in comparing Schlieren 
photographs with pressure distribution measurements. 

Errors in the pressure distribution measurements can 
arise from faulty orifices, lack of two-dimensional flow, 
etc. The results presented were checked far freedom 
from systematic errors by a simple procedure. The 
orifices were purposely distributed asymmetrically with 
respect to both centerlines of the airfoil. Rotating the 
airfoil through an angle of 180° moves each orifice into 
a new position. If the pressure distribution in both 
positions checks within the experimental scatter, the 
systematic errors are negligible. This was found to be 
the case for the measurements presented here. 


RESULTS 


*‘Laminar”’ and ‘‘Turbulent” Shock Waves 


Figs. 7 to 12 show a series of pressure distribution 
and Schlieren photographs for the airfoil with both lam- 
inar and turbulent boundary layer. It is evident that 
marked difference in the shock-wave pattern and pressure 
distribution exists between the laminar and the turbulent 


cases.* 

The main observations are these: The “laminar” 
shock waves are inclined upstream. No sudden in- 
crease in the surface pressures is observed. The ‘‘tur- 


bulent’’ shock waves, on the other hand, are nearly 
normal or inclined downstream. The pressure distri- 
butions show the familiar sudden pressure increase. 
In the laminar case two or more shock waves appear; 
in the turbulent case only one shock wave is present, in 


general. 


Shock ‘‘Reflections’’ from the Boundary Layer 


Close examination of a Schlieren photograph of the 
shock wave for the case of laminar boundary layer shows 
the following (Fig. 13): The shock wave is followed 
by a fan-shaped expansion zone. The boundary layer 
thickens ahead of the shock wave. At the point of con- 
tact with the shock wave and expansion zone, the edge 


* Similar ‘results have been found independently by Dr. J. 
Ackeret, at Zurich, Switzerland. Ackeret’s results apparently 
have not yet been published. Dr. H. S. Tsien who visited the 
Eidgenossische Technische Hochschule in 1945 informed the 
author of Dr. Ackeret’s investigations. 
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Fic. 7a. Pressure distribution 3-in. by 12 per cent circular arc airfoil, R = 8.33 X 105 (laminar). 


Fic. 7b. Flow at MM. = 0.795, R = 8.33 X 10° (laminar) 


Fic. 7c. 


Flow at 1. = 0.795, R = 8.33 X 105 (laminar); 





; knife edge normal to mean flow directon. 





knife edge parallel to mean flow direction. 
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Fic. 8a. Pressure distribution 3-in. by 12 per cent circular arc airfoil, R = 8.45 X 105 (laminar). 





Fic. 8b. Flow at MM, = 0.843, R = 8.45 X 105 (laminar); knife edge normal to mean flow direction. 


of the boundary layer is seen to make a definite bend 
toward the wall. 

The flow field existing near the base of the shock wave 
is shown schematically in Figs. 14 and 15. The flow 
passes the shock wave first and the pressure increases 
from P; to P,. The Mach Number decreases from M, 
to M2. Since the shock wave is oblique, the direction 
of flow changes by an amount A, = 8 — &. These 
changes are discontinuous—i.e., extremely rapid. 

In the subsequent expansion zone the flow is again 
accelerated, but in a continuous fashion, corresponding 
to the well-known Prandtl-Meyer flow. The pressure 


decreases from P, to P;; the Mach Number increases 
from M, to M3. The flow is turned farther in the same 
direction by an amount Ao’ = tH — ds. 

Quantitative results can easily be obtained from the 
Schlieren photographs using the results of the well- 
known theory of oblique shock waves and of Prandtl- 
Meyer flow (i.e., isentropic flow around a corner). 
The procedure is as follows: 

(1) Ms; is obtained from the angle @; between the 
direction of flow (indicated by the visible edge of the 
boundary layer) and the downstream edge of the ex- 
pansion zone. 63 is the Mach angle and, hence, 
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Fic. 8c. Flow at M. = 0.843, R = 8.45 X 10° (laminar); knife edge parallel to mean flow direction. 
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abel 


Meyer theory. 


My, is obtained from the angle of spread, w, of 
the expansion zone and M3, by using the Prandtl- 
(See, for example, reference 5.) 


wave. 


Mach Number, M2, behind the wave determine MJ, 
the Mach Number on the upstream side of the shock 
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.9b. Flow at M.. = 0.895, R = 8.77 X 10° (laminar); knife edge normal to mean flow direction. 


Fic. 9c. Flow at M. = 0.895, R = 8.77 X 105 (laminar); knife edge parallel to mean flow direction. 
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Fic. 10a. Pressure distribution 3-in. by 12 per cent circular arc airfoil + !/1»- by 3-in. flat plate, R = 1.67 X 10° (turbulent). 





Fic. 10b. Flow at M. = 0.795, R = 1.67 X 10 (turbulent); knife edge normal to mean flow direction. 
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Fic. 10c. Flow at M, = 0.795, R = 1.69 X 10° (turbulent); knife edge parallel to mean flow direction. 
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(4) From the stagnation pressure p) ahead of the 
shock and the Mach Number, the stagnation pressure 
behind the shock p)’ may be computed. 

Table 1 presents the results of a typical evaluation of 
this kind. 

It should be remembered that ordinary Prandtl- 
Meyer flow past a convex surface is an exact solution 
only for a supersonic zone of infinite extent. However, 
it can be expected to give a good approximation to a 
turning supersonic flow only when the radius of curva- 
ture is small compared with the extent of the supersonic 


It is evident that: 


(1) The flow behind the shock wave—expansion 
zone combination is supersonic. 
process can occur. 
subsonic speeds are finally obtained. Recompression 
occurs in two ways. 
shock wave is seen to increase slightly. The Schlieren 
photograph shows a faint triangular compression zone. 
The entropy increases in the passage through the shock 
wave, and thus the Mach Number corresponding to 
the same local pressure is less behind the shock. 


Hence another similar 
The process will be repeated until 


The pressure upstream of the 
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Fic. llb. Flow at VW, 


= 0.843, R = 1.69 X 10® (turbulent); knife edge normal to mean flow direction. 
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Fic. lle. Flow at /.. = 0.843, R = 1.69 X 10° (turbulent); knife edge parallel to mean flow direction. 
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Fic. 12a. Pressure distribution 3-in. by 12 per cent circular arc airfoil + '/,,- by 3-in. flat plate, R = 1.75 XK 10° (turbulent). 


(2) The pressures just upstream and just down- 
stream of the shock-expansion combination do not differ 
greatly. Hence, this type of shock wave is not easily dis- 
cerntble from surface pressure distribution measurements. 
This fact, which was first found by an analysis of 
Schlieren photographs, was completely borne out by 
the pressure distribution measurements made at a later 
date (see Fig. 8). The behavior of the pressure is 
sketched in Fig. 15. 

(3) The boundary layer does not separate at the 
shock-wave position. On the contrary, the flow turns 
toward the surface. 


The behavior of the shock wave near the surface can 
be summarized in the statement that the shock wave in- 
teracts with the boundary layer in a manner similar to 
shock interaction with a free jet surface. In other 
words the boundary condition satisfied here is more 
nearly the condition of constant pressure than that of 
constant direction. ‘This fact is believed to be of funda- 
mental importance for the understanding of many trans- 
onic phenomena. 

This result is by no means surprising. The rigorous 
boundary condition in any case is the condition of zero 
slip at the solid wall. This condition can be satisfied 
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only by a viscous fluid. The classical boundary layer 
theory leads to the result that the flow about a body 
can be separated into a boundary layer where viscous 
effects are appreciable and potential flow outside of this 
layer. The boundary condition of the viscous flow in 
the boundary layer is the condition of zero slip. The 
condition to be satisfied by the potential flow outside of 
thé boundary layer is that of zero normal velocity near 


the wall—i.e., at a distance equal to the displacement 


FLOW 


Flow at /.. = 0.895, R = 1.75 X 108 (turbulent); knife eige normal to mean flow direction. 


Flow at /. = 0.895, R = 1:75 X 10® (turbulent); knife edge parallel to mean flow direction. 


thickness of the layer. Hence, according to the classical 
boundary layer theory, potential flow is similar to flow 
past a solid body of a shape nearly identical with the 
real body. 

This result is based on the assumption that all changes 
of fluid properties in the direction of flow are small 
compared with the changes normal to the flow direc- 
tion (i.e., the boundary layer thickness 6 is assumed to 


be the smallest characteristic length). This condition 
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Fic. 13. Flow at M. = 0.843, R = 8.45 X 105, showing the flow expansion and the ‘‘corner’’ of the boundary layer; knife edge 
normal to mean flow direction. 


can fail near the base of a shock wave. The changes 
of velocity in the x direction—i.e., in this case parallel 
to the mean flow—are often of the same order or larger 
than the changes in the y direction—i.e., across the 
boundary layer. A new length, 7, the thickness of the 
shock wave, is important. Since r < 6, all results of the 
classical boundary layer theory become questionable in 
the neighborhood of a shock wave. This includes, of 
course, the boundary condition to be satisfied by the 
flow outside the viscous zone. If the changes in the di- 
rection of flow become large compared to the changes 
across the layer, the appropriate boundary condition 
for the outer flow appears to approach the condition pre- 
vatling near a free surface rather than near a solid wall. 
Although the inclination of the “laminar’’ shock wave 
upstream may appear surprising at first, it can easily 
be explained. Essentially, it is due to the finite extent 
of the supersonic zone. In terms of Mach lines—i.e., 
expansion and compression wavelets—the existence 
of a finite supersonic region means that the effects of 
both families of Mach lines must be considered. In 
general, for flow past a convex surface, the expansion 
wavelets will be inclined downstream and the compres- 
sion wavelets inclined upstream. Thus, if a shock wave 








Fic. 14. Nomenclature of shock wave-expansion zone. 


is formed as an envelope of compression wavelets, it 
must be inclined upstream. 


The Shock Wave in the ‘“‘Turbulent’’ Case 


Close examination of the shock-wave pattern for the 
turbulent case (Fig. 11) reveals a completely different 
result from that of the laminar case. The base of the 
shock wave is wide; thus, compression near the edge 
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Fic. 15. Pressure distribution in the neighborhood of a shock 
wave-expansion zone. 
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of the boundary layer is more nearly continuous. Ac- 
cording to the previous discussion, in these cases the 
flow should behave more in the conventional manner— 
i.e., corresponding to the classical boundary layer rela- 
tion. This is borne out by the pressure distributions 
(Figs. 10 to 12). For this case, a marked increase in 
pressure is evident at the shock-wave position. 


Separation and Shock Wave 


The argument presented in the preceding paragraphs 
concerning the application of results of the classical 
boundary layer theory also applies to considerations of 
boundary layer separation. The oft-repeated argu- 
ment, that the strong pressure gradients associated 
with a shock wave must cause separation, should be ap- 
plied with care. The investigations presented here 
show that separation does not always occur near the 
shock-wave position. This is, of course, evident in 
the case of the shock-wave-expansion zone combination. 
But even in the case of a strong single wave, such as the 
turbulent shock of Figs. 10 to 12, separation does not 
occur at the shock position. In fact, comparison of 
Figs. 11 and 12 shows the separation at approximately 
the same chordwise position, in spite of an appreciable 
shift in the location of the shock wave. A pronounced 
separation is clearly observed in the case of laminar flow 
at M = 0.895. Here, the boundary layer separates 


ahead of the shock wave. A schematic drawing of this 


condition is presented in Fig. 16. 


CONCLUSIONS 


(1) The boundary layer has a pronounced influence 
upon the shock-wave pattern and pressure distribution 
in transonic flow. A change of the boundary layer 
from the laminar to the turbulent régime results in a 
marked change in shock-wave pattern and surface 
pressures at the same free stream Mach Number. 

The importance of this finding for the interpretation 
of wind-tunnel results is evident. Effects of this type 
should also have a considerable influence on the pres- 
sure distribution on turbine and compressor blades. 

(2) Shock waves can interact with the boundary 
layer in a manner similar to the reflection of a shock 
wave from a free surface. The shock wave in this case 











Fic. 16. Boundary layer separation ahead of a shock wave. 


is followed by an expansion zone; the surface pres- 
sure does not show a marked change near the shock- 
wave base. The pressure gradient normal to the bound- 
ary layer is of the same order as the pressure gradient 
parallel to the boundary layer. 

This effect is in striking contrast to the behavior of 
boundary layers in subsonic flow and is believed to be 
important for an understanding of supersonic flow near 
solid surfaces. 

(3) It is found that separation does not necessarily 
occur near the base of a shock wave. 

Arguments that postulate the separation due to a 
shock wave are based upon the classical boundary layer 
theory, which does not apply in the immediate neigh- 
bornood of the shock-wave base. 
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Letters to the Editor 


Dear Sir: 

E. V. Laitone’s Letter, to the Editor, published in the August, 
1946, issue of the JouRNAL (p. 404), is concerned with a subject 
upon which the present writer has already commented elsewhere.! 
It was pointed out that there are several variants of the linear- 
perturbation theory, and by different writers, all of which pro- 
vide the same results in two-dimensional cases (such as the cor- 


rection of airfoil characteristics for compressibility effects), but 
that these variants lead to inconsistent and incompatible results 
when applied to bodies of revolution in axially symmetric flow. 
Unfortunately, this somewhat confused subject was rendered a 
bit more obscure by the writer’s note in the Quarterly, which in- 
cluded the statement that all of these results in the axially sym- 
metric case are equally valid and cannot be distinguished between 
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in a linearized theory. This statement is erroneous, since further 
investigation has disclosed that the liberties taken with the 
boundary condition in the plane case (specifically, fitting the 
boundary values of the transverse velocity at an approximate 
location) are not permissible in the three-dimensional case, and 
only one of the variants is innocent of this error. Thus the veloc- 
ity-correction formula loses its ambiguity, and its only correct 
expression is as follows: 

(u — U)/U = (1/6*) F(6n) (1) 
where u is the surface velocity, U is the velocity of the stream, 
8? = 1 — M? (M being the stream Mach Number), 1 is the 
parameter defining the slenderness of the body, such as the ratio 
caliber/length, and F(m) is the value of (u — U)/U at the corre- 
sponding point of a body of slenderness m in incompressible flow 
(8 = 1). This is the result given by B. Géthert.2, The writer 
has prepared a second note for publication in the Quarterly which 
will present these matters in more detail. 

Thus, Laitone is correct in his objections to the, use of the 
velocity-correction formula F(m)/8. Moreover, he is correct in 
his suggestion that the surface velocity is approximately inde- 
pendent of Mach Number for very slender bodies, because for 
typical bodies 


F(n) = n? log n 


and 
(1/8?) F(@n) = n? log (8n) = n* log n 


provided n < 8. But this approximation is not very good for 
bodies of practical interest, for which 1/10 < » < 1/3 and 0 < 
M S 0.8, say. It is recommended that Géthert’s formula, Eq. 
(1) be used; it results in velocity corrections as great as 40 per 
cent when applied to ellipsoids of revolution in this régime. 

It is of interest also that Sauer’s proof that the velocity is inde- 
pendent of 8 in the limiting case, which Laitone quotes, is really 
invalid, since it involves an improper integral that does not ac- 
tually exist. This has been pointed out by Dr. C. C. Lin, of 
Brown University. Lin has provided a valid proof of the theorem 
and has arrived at a formula which shows (u — U)/U to be pro- 
portional to nm? log (8m), for small n, in general. Previously this 
form was obtained by Géthert? and by the writer! only for the 
special case of ellipsoids of revolution. 
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Dear Sir: 
In a Letter to the Editor which appeared in the August, 1946, 
issue of the JouRNAL, Edmund V. Laitone pointed out that the 


“correction factor” 1/171 — M? is applicable only to two-di- 
mensional subsonic compressible flows and not to three-dimen- 
sional flows, as discovered by Ferrari! and also by Géthert.? 
However, it is not correct to state that the pressures acting on a 
slender body of revolution are independent of Mach Number in 
first approximation. A more careful analysis* shows that while 
the equivalent source-sink distribution for a slender body of 
revolution is independent of Mach Number in first approximation, 
the increase with Mach Number of the maximum (positive) 
axial disturbance velocity, or the maximum negative pressure 
coefficient on such a body, is approximately given by the follow- 
ing expression: 

Com _ y 4 log V1 — M? (1) 


Cpo 0.31 + log t/c 


where ¢t/c = thickness ratio of the body. (This relation, which 
was derived by considering a slender ellipsoid of revolution, agrees 
with the result obtained by Schmieden and Kawalki‘ by a strict 
application of Géthert’s method.) Thus, Mr. Laitone’s state- 
ment is correct only for the limiting case of a ‘“‘body”’ of zero 
thickness. 

Apparently, Sauer® and others are led to the conclusion that 
the pressures are independent of Mach Number because in their 
analysis the contribution of a source or sink element to the axial 
disturbance velocity, which is proportional to 


(x —é) f(de 
(x — 8? + (1 — Mr? 





is replaced by f(&)/(x — £)dé. Since the term (1 — M?)r? is dis- 
carded, the effect of Mach Number on the pressures is also dis- 
regarded. It is difficult to justify such a procedure, which 
amounts to neglecting the logarithmic term in Eq. (1) in compari- 
son with unity. In references 3 and 4 the solution of the linear- 
ized partial differential equation for the velocity potential is ac- 
cepted without further modification, except that in the final re- 
sult terms of the order of (¢/c)? log (t/c) are neglected in compari- 
son with terms of the order of t/c. 
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Jet Propulsion Applied to Helicopter Rotors 


R. H. MILLER* 
McDonnell Aircraft Corporation 


SUMMARY 

Some of the factors affecting jet-driven rotor performance are 
discussed and it is shown that the relatively high tip speeds re- 
quired in order to realize a reasonable jet efficiency result in a 
considerable reduction in the aerodynamic efficiency of the rotor. 
Although this aerodynamic efficiency may be improved by de- 
creasing the rotor solidity, a practical limit is reached, in the 
case of jet drives involving airflow through the bladés, beyond 
which the specific fuel consumption increases almost asymp- 
totically with decrease in solidity. However, the weight-saving 
possible with such drives, and their simplicity, may in certain 
cases offset their low aerodynamic and propulsive efficiencies. 
A rough comparative performance and weight estimate are made 
between conventionai and jet-driven helicopters in an attempt to 
evaluate their relative merits. 


INTRODUCTION 


— POSSIBLE ADVANTAGES of a jet-propelled rotor 
are (1) elimination of antitorque devices, and (2) 
weight-saving and mechanical simplicity. However, 
in order to obtain a reasonable jet propulsive efficiency, 
the rotor should be operated at relatively high tip 
speeds, which will result in a reduction in the aerody- 
namic efficiency of the rotor unless the solidity is ma- 
terially decreased. The minimum rotor solidity, in 
the case of jet drives involving air or gas flow through 
the blades, is in turn governed by the necessity 
for keeping the gas velocities in the blade, and hence 
the associated pressure’ drops, to a _ reasonable 
value. 

Several different jet drives for a helicopter rotor may 
be envisaged but it is believed that the problems com- 
mon to most of them may be investigated by consider- 
ing only the following two basic types: 

Type A.—Propulsion by means of a jet with com- 
pression provided entirely by virtue of the forward 
motion of the jet unit located at or near the blade tip, 
i.e., a ram-jet. 

Type B.—Propulsion by means of a jet at the blade 
tip with centrifugal compression in the hollow blade and 
additional compression provided by an engine or rotor- 
driven compressor. 

Obviously in the case of Type A, the jet efficiency is 
independent of the rotor solidity except as structural 
considerations limit the-rotor blade chord for a given 
jet unit weight. In the case of Type B, the jet effi- 
ciency is extremely sensitive to rotor solidity. 


Received March 14, 1946. ' 
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NOMENCLATURE 


A = disc area, sq.ft. 
C = blade chord 
c = specific fuel consumption (S.F.C.), lbs. per b.hp.-hr. 


Cp = specific heat at constant pressure = 0.24 
Cr = thrust coefficient = thrust/pAQ?R? 
J = mechanical equivalent of heat = 778 ft.lbs. per B.t.u. 


M, = weight flow of air, lbs. per sec. 
My = weight flow of fuel, Ibs. per sec. 


P = power, ft.lbs. per sec. 

R_ = blade radius, ft. 

T = temperature, °F. absolute 
V = tip speed, ft. per sec. 


W = gross weight, lbs. 
or work, ft.lbs. per sec. 
f = equivalent flat plate area—parasite drag 
or Fanning friction factor 
= acceleration of gravity, ft. per sec.? 


g 

k = ratio of specific heats = 1.4 

r = presstire ratio 

“u = jet velocity, ft. per sec. 

w = disc loading, lbs. per sq.ft. 

6 = profile drag coefficient of rotor blade 

o@ = rotor solidity—ratio of blade area to disc area 

Q = rotor speed, rad. per sec. 

p = density of air = 0.002378 lbs. sec.? per ft.‘ 

n = efficiency of jet drive 

7, = mechanical efficiency of rotor drive: ratio of shaft 
to brake horsepower 

n. = thermal efficiency of engine 

# = ratio of forward speed to tip speed 


SECTION I. AERODYNAMIC EFFICIENCY OF THE ROTOR 


In order to permit an evaluation of the effects of 
changes in tip speed and solidity on the aerodynamic 
efficiency of a rotor, the weight that may be lifted for 
a given shaft horsepower and the maximum range for 
a given ratio of fuel load to gross weight will be estab- 
lished as a function of solidity and tip speed. 

The well-known formula for the hovering helicopter, 
assuming constant profile drag, 5, gives for the torque 
coefficient 


Ce = 05/8 + Cr C,/2 


which may be expressed in terms of power and weight as 


<n (2) Q3R3 + y- 
W 8 \w 2p 
This has been plotted in Fig. 1, in terms of Ibs. per 
shaft hp., against o for various tip speeds and disc 
loadings. 


Additional information as to the vertical rate of climb 
for a given excess power, AP (in ft.lbs. per sec.) over 
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Fic.1. Effect of rotor solidity, o, tip speed, QR, and disc loading, 
w, on power required to hover. 
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Fic. 2. Ratio of forward speed to tip speed, », for maximum 
range for various thrust coefficients and solidities. 


that required for hovering as obtained from Fig. 1, 
may be readily determined, for a constant-speed rotor, 
from the relationship 


Vertical rate of climb = 





AP on + _—) 
W \AP/W + Vw/2p 


in ft. per sec. which may be derived from the standard 
formulas for vertical climb. 

The range at any forward velocity, .QR, may be 
closely approximated from the formula 


Range = 
375u ACr (n/c) 





(05/8)(1 + mu2) + (Cp*/2u) + (f/A)(u8/2) in miles 


m may be taken as 4.6 and f/A, the ‘‘drag coefficient” 
based on disc area for the helicopter less rotor, as 0.010. 


AC, is the increment of thrust coefficient represented 
by the fuel. C, should be based on the average gross 
weight. 

The value of » for maximum range may be established 
by solving for uw in the formula 

o5 _ (Cr*/u) — (wif/A) 

8 nu? — 1 
the solution of which has been plotted in Fig. 2 for a 
drag coefficient f/A of 0.010 and a value of 6 of 0.012. 

The maximum range has been determined for various 
values of tip speed, disc loading, and solidity, and 
plotted against solidity in Fig. 3. These curves, and 
those of Fig. 1, have been stopped at a minimum rotor 
solidity corresponding to a blade angle of attack of 10°, 
Beyond this angle of attack, which is believed to be 
the maximum desirable operating angle of attack, the 
drag coefficient increases rapidly and the assumption 
of a constant drag coefficient, 6 = 0.012, no longer 
holds. Fig. 3 is based on a fuel load of 7.5 per cent 
of the gross weight, but for any other fuel load it may 
be assumed, as a close approximation, that the range 
increases proportionately to the increase in fuel. It 
is evident from Figs. 1 and 3 that operation at higher 
tip speeds entails a considerable loss in both cruising 
and hovering performance unless the solidity is de- 
creased accordingly. 

However, the reduced power-plant and transmission 
weights and the elimination of the antitorque rotor 
result in a considerably lighter helicopter when a jet 
rotor drive is used. It would therefore appear possible, 
by increasing the fuel load and shaft horsepower, to 
equal or exceed the performance of a conventional 
helicopter without sacrificing the load-carrying ca- 
pacity. This will be investigated in the next sections. 

Losses due to compressibility effects, which may be 
appreciable at the higher tip speeds and lower solidities, 
have been neglected in the preparation of the curves 
of Figs. ] and 3. 





Section II, Jet EFFICIENCIES 


Ideal Cycle Efficiencies 


For any jet operating at a pressure ratio r (defined 
as the ratio of combustion chamber pressure to atmos- 






































. ” =.60 — wae 
w . 
8 c =-60 --- wal s 
= 20 i b = 
; Whang Sai: meee etd ein | 
z «. ~*~ | ——_[ } ar-lese 
100) — = _ a=. 
: le i ee 
~~ -2R=900 
° -02 04 -06 08 10 +2 14 


cs 
Fic. 3. Effect of rotor solidity, 7, tip speed, QR, and disc 


loading, w, on maximum range. Specific fuel consumption, ¢ 
assumed to be 0.6. 
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pheric pressure) the work of compression per pound of 
air is 


W. = JtpTo(r*-?/* — 1) 


where 
J = mechanical equivalent of heat—778 ft.lbs. per 
B.t.u. 
Cp = specific heat at constant pressure = 0.24 
k = c,/c, = 1.4 for air at standard temperature 


whence (k — 1)/k = 0.286, which will be used hereafter. 
Ty) = ambient air temperature = 520°F. absolute 


This work of compression may be supplied by any 
means and may be assumed independent of the rotor 
speed. The work done in accelerating one pound of 
air to the rotor tip speed, V, is 

W, = V?/2g 

The total work available from the jet per pound of 

charge air is 
W, = uV/g 
where u is the jet velocity. If the velocity in the 


combustion chamber is assumed to be negligible in 
comparison with V, then 








“= V—2¢FcT2 [(1/r-*) — 1] 


where 7% is the combustion chamber temperature. 

For combustion at constant pressure, the tempera- 
ture of the gases in the combustion chamber, T>, is 
given closely by 


My/Ma4 = ¢(T2 — T1)/18,000 
where 


Mr weight flow rate of the fuel 
Ma = weight flow rate of the air 


Il 


18,000 = heat content of the fuel, including an 
overall combustion chamber efficiency 
of 90 per cent 

Ti = temperature of air entering the combustion 


chamber = 7or°-786 


The total work available in the fuel is, assuming 
20,000 B.t.u. per Ib., 


W, = JMp(20,000) 
= Jc,M4(T2 — T;)(20,000/18,000) 


and the total work available per pound of charge air 
is then : 


W; = Jc,(1/0.90)(T2 — T:) 


When compression is provided entirely by the rotor, by a combination of centrifugal compression in the blades 
and a rotor-driven compressor, then the overall efficiency of the drive is given by 


n = (W; — W, — W.)/Wr 





V2 


cal 1 
24 —2gJc,T2 (Ga- 1) - > — Je,To(r-** — 1) 





am 208(T: — 7) 


This relationship has been plotted against r for 
various values of 72 in Figs. 4a, 4b, and 4c for rotor 
tip speeds of 600, 900, and 1,120 ft. per sec. If the 
work of compression is supplied only by ram or cen- 
trifugal compression in the rotor blade, the’ work of 
compression, W, and W, are 


W. + We = V*/g 


The pressure ratios corresponding to pure ram 
compression are indicated in Fig. 4 and these points 
represent the efficiencies of jets of Type A. 

The favorable pressure distribution on the forward 
end of a ram-jet resulting from the external airflow 
induced by the jet may increase these efficiencies by 
a considerable amount but no attempt is made here 
to evaluate these effects. 

If additional compression is provided by means of an 
engine-driven compressor, rather than by a rotor- 





driven compressor, the overall efficiency is given by the 
expression 











JtpTo 
Ne 


208(72 ~ T1) + (1719-286 — 1) 


where V?/g is the work of acceleration and compression 
in the blade and 7 is the pressure ratio across the 
engine-driven compressor. If r2 is the compression 
due to centrifugal forces in the blade, then r = fori. 
Ne is the overall thermodynamic efficiency of the engine 
used to drive the compressor. 

Figs. 5a, 5b, and 5c show the efficiency plotted 


against pressure ratio fer this combination, assuming 
a thermodynamic efficiency for the engine of 25 per cent 
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Fic. 4. Theoretical efficiency of Type B jets, with additional 
compression provided by a rotor-driven compressor, for various 
combustion chamber temperatures, 72, and pressure ratios, r. 
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_ Fic. 5. Theoretical efficiency of Type B jets, with additional 
compression provided by an engine-driven compressor, for 
various combustion chamber temperatures, 72, and pressure 
ratios. 7. 
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corresponding to a specific fuel consumption of 0.6 
Ib. per b.hp.-hr. 

Figs. 4 and 5 cover all possible combinations of the 
two types.of jet drives, A and B, listed above. 

All of the above efficiencies are theoretical ideal 
efficiencies, i.e., mo compressor, turbine, or nozzle 
losses have been considered, although a combustion 
efficiency of 90 per cent has been included in order to 
keep combustion and cycle losses separate. 

In addition, the following assumptions, which it is 
believed have little effect on the final results, have been 
made: 

(1) The velocity of the air leaving the compressor 
and in the rotor blades is small compared to the velocity 
of the jet and the velocity of the blade tip. 

(2) The value of c, is constant and equal to 0.24. 
This neglects the variation of c, with temperature and 
gas content. The same assumption has been made 
with regard to k. 

(3) Combustion occurs at constant pressure and 
with a negligible change in combustion chamber 
velocity. 

These assumptions, and further assumption of zero 
nozzle losses, have little effect on the final results. 

However, neglecting the pressure losses through the 
blades and combustion chamber has a considerable 
effect on the absolute value of the overall efficiencies 
shown in Figs. 4 and 5, although the relative values 
will not be appreciably affected. Provided the velocity 
of flow through the blade is kept below 200 ft. per sec., 
the blade and combustion chamber losses are negligible, 
but in order to maintain this velocity it may be neces- 
sary to provide a rotor with an excessively high solidity. 
The effect of these major losses on the specific fuel 
consumption will be investigated next. 


Effect of Losses on Specific Fuel Consumpton 


The main losses in efficiency associated with a jet 
rotor drive for a helicopter arise as the result of: 

(a) The necessity for pumping a _ considerable 
amount of air through the blades if a reasonable shaft 
horsepower is to be realized. 

(b) The necessity for keeping the dimensions of the 
combustion chamber located at the blade tip as small 
as possible in order to minimize the external drag. 
This increases the velocity of the gases in the combus- 
tion chamber and, consequently, the pressure drop 
across it. 

Loss (a) does not occur in the case of a ram-jet, but 
loss (b) may be more severe because of the low pressure 
ratio available. 

Since an increase in solidity for the purpose of re- 
ducing these losses to a negligible value results, beyond 
a certain point, in a considerable decrease in the aero- 
dynamic efficiency of the rotor, the limits between 
excessive losses in drive efficiency and excessive losses 
in aerodynamic efficiency must be established. 
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In order to estimate the lowest practical solidity for 
a given horsepower and rotor diameter, the losses in 
the blade have been estimated assuming a Fanning 
friction factor of 0.01, and the pressure drop through 
the combustion chamber has been assumed to be four 
times the dynamic pressure at the chamber inlet, 
with an inlet velocity equal to the velocity in the 
blade. 

In the case of jet B, an overall compressor efficiency 
of 75 per cent has been assumed, which is intended 
to include a compressor efficiency of 85 per cent, the 
maximum that may be realized with an axial flow com- 
pressor, and pumping losses of 10 per cent. These 
pumping losses arise as a result of the pressure drops 
associated with the airflow to the blades through the 
fuselage ducts and past the flap and lag hinges. 

In estimating the pressure drop through the blades, 
the density corresponding to the average temperature 
and pressure in the blade for the conditions being 
investigated has been used. 

In the case of the ram-jet (Type A), a combustion 
chamber velocity of one-quarter of the forward velocity 
and an entrance and diffuser loss of 5 per cent of the 
available ram have been assumed. Entrance losses 
have been neglected with jets of Type B because of the 
ram available due to the motion of the helicopter. 

Taking these losses into account, the specific fuel 
consumption for jets of Type B with a rotor-driven 
compressor and with an engine-driven compressor 
and jets of Type A has been plotted against the rotor 
solidity, o, in Fig. 6 for a tip speed of 600 ft. per sec. 
In addition, the fuel consumption for jets of Type B 
with engine drive and a ram-jet has been plotted for a 
tip speed of 900 ft. per sec. 

A 40-ft. diameter, two-bladed rotor has been assumed 
with an available net shaft horsepower of 150 hp. It 
is at once clear from an examination of Fig. 6 that a 
definite minimum rotor solidity exists for any type of 
jet drive using the blades as air passages. Since an 
increase in pressure drop results in a decrease in shaft 
horsepower, which in turn requires an increase in airflow 
resulting in a greater pressure drop, it is evident that, 
for a given shaft horsepower, when the critical solidity 
is reached the S.F.C. increases almost asymptotically. 
Operation at any point other than the flat portion of 
the curve would be, in practice, impossible. In the 
case of the jets of Type B with a rotor-driven com- 
pressor, the minimum solidity is far too high for aero- 
dynamic efficiency because of the necessity for supply- 
ing by means of the relatively inefficient rotor not 
only the 150 shaft hp, but also the power to drive the 
compressor. Such a type of rotor would therefore 
appear to be of little practical use for helicopters. 

While the curves of Fig. 6 have, of necessity, been 
prepared for a specific rotor diameter and horsepower, 
the effect of various changes may be roughly estimated 
from the following considerations. 
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Fic. 6. Effect of rotor solidity, ¢, on specific fuel consumption 
for various types of jet drives. Combustion chamber temperature 
in all cases is 2,000°F. abs. 


The available shaft horsepower is a direct function 
of the mass flow, M, of air through the blade. The 
pressure drop, Ap, through the blade is approximately 
a function of the square of the velocity, V, which, for 
given operating conditions, is a function of the mass 
flow divided by the cross-sectional area of the blade. 
The cross-sectional area of the blade is a function of 
the square of the blade chord, c. 

Therefore, 


Ap = f(V?) = F(M/c?)? 


whence for a given Ap, which is the major factor in 
determining the minimum rotor solidity, the maximum 
horsepower that may be supplied to the rotor varies 
as the square of the blade chord, or solidity. 

Thus, with a jet of Type B, r = 2, and an engine- 
driven compressor, if the solidity were increased from 
the minimum of approximately 0.04 to 0.06, the 
maximum permissible rotor horsepower could be 
increased from 150 hp. to 150 (0.06/0.04)? or 338 hp. 
without any appreciable increase in §.F.C. 

Also, if a single-rotor blade were used instead of two 
blades and the solidity were kept at 0.04, the cross- 
sectional area would double and the maximum per- 
missible horsepower would also double to 300 hp. 
Actually, the increase would be slightly greater because 
of the reduction in wetted area, as discussed below. 
Conversely, if the horsepower were held at 150 hp., 
the solidity could be reduced to 0.04(1/+/2) or 0.028. 

The pressure drop has been given above as a function 
of the chord to the fourth power. Actually, for the 
blade 


perimeter 





Ap = VopV7f X — X blade length 


cross-sectional area 


where f is the Fanning friction factor. The pressure 
drop along the blade, therefore, varies inversely as 
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TABLE 1 
Comparative Performance of Jet-Driven and Conventional Helicopters 
; Vertical 
Available Rate of 
Gross : Shaft Number Maximum Climb— 
Weight, Disc we Hp. at of Tip Fuel— Range— Ft. per Min. 
Lbs. Loading Solidity Hovering Blades Speed Gal. Miles at Sea Level Notes 
1 2,500 2 0. 06_ 140 3 450 30 192 0 Conventional helicopter 
2 2,250 1.8 0.045 150 2 600 68 127 0 Jet Type B engine-driven 
compressor 
3 2,330 1.9 0.032 150 1 600 80 192 370 Jet Type B engine-driven 
compressor 
4 5,000 2.55 0.06 350 3 450 75 275 800 Conventional helicopter 
5 5,450 2.80 0.06 505 2 600 285 275 1,710 Jet Type B engine-driven 
; - compressor 
6 5,150 - 2.62 0.043 505 1 600 240 275 2,480 Jet Type B engine-driven 
compressor 
7 1,060 3.00 0.020 100 2 900 30 80 1,000 Ram-jet-driven rotor 





¢ rather than c‘, if the variation of f with Reynold’s 
Number is neglected. However, the combustion 
chamber pressure drop, which varies more closely as 
c*, represents approximately half the total pressure 
drop between the blade root and the jet. Assuming 
Ap to be a function of c‘ is therefore slightly conserva- 
tive for increasing solidities and slightly nonconserva- 
tive for decreasing solidities. 

If the rotor diameter were increased and the solidity 
remained constant, the only pressure drop effected 
would be that due to friction in the blade which varies 
with c®. However, there is an increased pressure drop 
due to the increased blade length resulting in a net 
variation proportional to c‘, as before. Thus, if the 
rotor diameter were increased to 50 ft. and the solidity 
remained the same, the permissible shaft horsepower 
at the minimum solidity of 0.04 for the example taken 
above would become 150(50/40)? or 224 hp. If the 
solidity were increased to 0.06,.the permissible horse- 
power would become 224(0.06/0.04)? or 505 hp. The 
power required from the engine driving the compressor 
would increase directly as the increase in rotor power 
to 35(505/150) or 120 b.hp. 


SecTIon III. COMPARATIVE PERFORMANCE OF JET 
AND CONVENTIONAL HELICOPTERS 


In order to obtain approximate comparative per- 
formance estimates for a helicopter with a jet-propelled 
rotor and one with a conventional rotor drive, the 
following weight estimates will be used. 


(a) Conventional helicopter, 40-ft., three-bladed rotor, 
180-b.hp. engine: 

(1) Gross weight = 2,500 Ibs. 

(2) Fuel weight, including tanks = 200 lbs. 
(30 gal.). 

(3) Weight of complete power plant, less anti- 
torque device = 600 Ibs. 

(4) Weight of tail rotor, including supporting 
structure and transmission = 250 lbs. 

(5) Rotor group = 350 lbs. 


(6) Weight less items 2, 3, 4, and 5 = 1,100 Ibs. 


(b) Equivalent jet-propelled helicopter, with Type 
B jet, 35-b.hp. engine driving the compressor (see Fig. 
6)—using item (a6) above as a base: 


(1) Weight of power plant, compressor, gearing 
and ducts = 250 Ibs. 

(2) Rotor group = 450 Ibs. 

(3) Gross weight, less fuel: items (a6) plus 
(b1) and (b2) = 1,800 Ibs. 


(c) Conventional helicopter, 50-ft. rotor, 450-b.hp. 
engine: 


(1) Gross weight = 5,000 Ibs. 

(2) Fuel weight, including tanks = 500 Ibs. (75 
gal.). 

(3) Power plant, as in (a3) above = 1,250 lbs. 

(4) Tail rotor, as in (a4) above = 500 Ibs. 

(5) Rotor group = 550 Ibs. 

(6) Weight less items 2, 3, 4, and 5 = 2,200 Ibs. 


(d) Equivalent jet-propelled helicopter, with Type B 
jet, 125-b.hp. engine driving the compressor—using item 
(c6) above as a base: 


(1) Weight of power plant, as in (bl) = 700 lbs. 

(2) Rotor group = 650 lbs. 

(3) Gross weight, less fuel: items (c6) plus (d1) 
and (d2) = 3,550 Ibs. 


(e) Ram-jet-propelled helicopter—this will be esti- 
mated in a different manner because of the extreme 
simplicity of construction possible with this type of 
jet. The following weights are assumed: 


MT EAC MA liale ced wheel hleclaldd wake cts 150 Ibs 
Ns Adds ict. dd bb teehee dito ceded ons 200 Ibs. 
a iano ak casero epeieme 50 Ibs. 
Wiectricel Caetinsitemt. ...ccc ccc cotc vec ccces 40 Ibs. 
Comtrels. ...... 0. Fawitiehpiauana cco ua haul 4 ecocetere-een 20 Ibs. 
Ge ao 8k ch Oiiin ih oid he lds CeK4 Ses 400 Ibs. 
PRI, 5 os bo izbie nls oss eco ok rege Os 200 Ibs 
EE Sap. o occ 6500s veceeveesee 1,060 Ibs 
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The only limits on the minimum solidity of such a 
rotor are structural, since no air is pumped through the 
blades. A minimum practical solidity of 0.02 will be 
assumed. 

Table 1 has been prepared on the basis of the fore- 
going weight estimates, the aerodynamic performance 
data of Figs. 1 and 3, and the specific fuel consumption 
of Fig. 6. The effect of changes in disc loading and 
fuel load has been estimated by interpolation and 
extrapolation of Figs. 1 and 3. The effect of changes 
in solidity, number of blades, and diameter from that 
assumed in Fig. 6 has been estimated as desctibed in 
Section IT. 

It is evident from Table 1 that except for the ram- 
jet (7) it will probably be possible in every case to 
exceed the performance of a conventional helicopter, 
provided that a single-blade rotor could be built and 
made to operate satisfactorily. If a two-bladed rotor 
is used, the performance will be still better, except for 
the case of the light helicopter with a jet rotor of 
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Type B, which has less range than the conventional 
helicopter. 

Although the fuel load of the jet helicopter is in all 
cases considerably greater than that of the conventional 
helicopter, this does not necessarily represent an 
increase in operating costs because of the much lower 
grades of fuel that may be used. 

In estimating the performance of the conventional 
helicopters, it has been assumed that, after deducting 
cooling fan, drive, and antitorque losses, the net shaft 
horsepower remaining is 77'/2 per cent of the available 
engine brake horsepower in hovering and 85 per cent 
in cruising. 

Several other combinations of solidity, fuel load, 
shaft power, and tip speed other than those investi- 
gated and presented in Table 1 are possible and may 
result in even better jet rotor performance, although 
at the expense of an increase in gross weight. How- 
ever, it is believed that the general trend is clearly 
indicated by the estimates of Table 1. 





Letter to 


Dear Sir: : 

From the Schlieren photographs and discussion given by Stack 
in the Eighth Wright Brothers Lecture,' the writer developed 
the boundary layer-shock interaction theory explained in ref- 
erence 2. 

The flow patterns shown in Fig. 1 may be constructed as shown 
in reference 2 from the following considerations: 

(a) The Mach waves or infinitesimal oblique shocks are 
inclined to the local velocity vector at the angle 


Oy sin~!(1/Myz) (1) 


Laminar Boundary Layer 


the Editor 


Turbulent Boundary Layer 


(b) The local supersonic isentropic flow compression (H, 
He) produced in the contraction of the stream tube above the 
thickening boundary layer is given by 





b1/He c+  —— P 
— s& = 06 
pi/H= a ld 
24 
Y— Ms — yMit + yMi 
bi Qs? (2) 
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Fic. 1. Effect of increasing boundary layer turbulence or free-stream Reynolds Number. 
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X/¢ 
Fic. 2. NACA 4412 airfoil, M = 0.77,a = —0.25°. 


where Heo is the free stream total pressure and p; and Mj, are 
the static pressure and Mach Number at the point on the airfoil 
surface where the sudden thickening of the boundary layer starts 
(@g = 0). The angle @g in radians is the relative deflection of 
the local velocity vector immediately above the boundary layer 
as the boundary layer thickness increases. Eq. (2) may be 





written as 

¥+1 1 1 
‘MEP ee oho = we * lh...) 
. ym \p p} 





tv Mi —] 
(3) 


Fig. 1 indicates that the sudden -boundary layer thickening 
produces isentropic flow compression and decreases the super- 
sonic flow velocity until the local Mach Number (1/,) approaches 
unity or the boundary layer cannot separate any further. The 
final flow compression to a subsonic flow is then accomplished 
by means of a small normal shock. As the amount of boundary 
layer separation which the given flow can produce is reduced, 
the final normal shock increases in magnitude as indicated in 
Fig. 1. 

Although the flow above the boundary layer (dotted line) has 
become subsonic behind the normal shock, Fig. 1 indicates that 
the flow inside the thickened boundary layer (solid line indicating 
pressures on airfoil surface) may still be supersonic. This is not 
so, howeyer, as may be seen by an analysis in Fig. 2 of the total 
pressure measurements inside the boundary layer. The total 
pressure loss given at x/c = 1 was taken from reference 3 as the 
measurement near the airfoil trailing edge. The total pressure 
loss inside the boundary layer was then faired in accordance 
with the derivation in reference 2. The experimental total 
pressure loss through the normal shock was in exact agreement 
with that predicted by theory so the variation of total 
pressure was taken as the sudden discontinuity predicted by 
theory. 

Fig. 2 shows that the actual local Mach Number inside the 
boundary layer has become subsonic due to the total pressure 
loss from viscosity effects. It is seen that throughout the 
isentropic supersonic flow compression the static pressure ), 
is continuous as in ordinary boundary layer theory but only up 
to the normal shock. Downstream of the normal shock there 
is a pressure discontinuity between the flow inside and outside 
of the boundary layer in agreement with the Schlieren photo- 
graphs of reference 1. However, the static pressure inside the 
boundary layer is less, not greater, than that of the flow outside 
the boundary layer. 

Consequently, it is not necessary to assume that only the 
sudden pressure gradient of a shock wave produces the abrupt 
boundary layer thickening. Actually, if the boundary layer 
Reynolds Number is such that an infinitesimal pressure increase 
produces a thickening of the boundary layer, then a Mach wave 
or infinitesimal oblique shock arising from any temporary dis- 
turbance can start the thickening of the boundary layer and the 
resulting isentropic flow compression that will then continue 
until it finally concludes in a small normal shock. If the bound- 
ary layer is sensitive to an increase in the pressure gradient, then 
the initial thickening will be abrypt and the initial Mach waves 
or infinitesimal oblique shocks will tend to converge and form a 
finite oblique shock. However, this will still be a relatively 
weak oblique shock so that Eqs. (1), (2), and (3) may still be 
used to compute the flow relations. 
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Newton and Fluid Mechanics 


J. C. HUNSAKER* 


National Advisory Committee for Aeronautics 


Pees CONCEPTION OF DYNAMICS armed his suc- 
cessors of the next three centuries with basic tools 
with which they erected the great structure of modern 
engineering science. : 

While Newton is claimed by mathematicians and as- 
tronomers as their own, engineers owe to him the very 
foundation of their art. Modern aerodynamics is com- 
pletely Newtonian in its development. Not only is it a 
consequence of Newton's laws of motion, but to this day 
it continues to utilize some of his original tactics to ob- 
tain solutions. 

For example, Newton, being unable to determine 
from first principles the resistance of a body moving 
through a real fluid, simplified his problem by separating 
it into three parts. He postulated, first, a special fric- 
tionless incompressible fiuid, then a viscous fluid, and 
finally, a compressible fluid. 

For the frictionless fluid, he deduced that there would 
be a resistance due to the impact of fluid particles, vary- 
ing as the fluid density, as the square of a linear dimen- 
sion of the body, and as the square of the velocity of 
motion. 

For the second fluid, having ‘‘a want of lubricity’’ or 
viscosity, he concluded that the frictional force must 
vary as the rate of shear of adjacent layers of fluid. 
This is a clear definition of the coefficient of viscosity 
and the basic characteristic of laminar flow. To this day 
we speak of tars and greases as non-Newtonian fluids, 
because they do not exhibit this linear relation between 
force and rate of shear. 

Finally, he speculated on the propagation of pressure 
pulses, like sound waves, in a compressible fluid and 
found the velocity of propagation to be a function of the 
density and elasticity of the fluid. His computation for 
the velocity of sound in air was a fair approximation 
which stood until La Place corrected it for adiabatic 
rather than isothermal conditions. : 

From the concept of the ideal frictionless fluid comes 
the classical hydromechanics of Bernoulli, Euler, 
d’Alembert, and Lagrange. . 

From the concept of the viscous fluid we have the 
great development by Stokes and Osborne Reynolds 
leading to Prandtl’s boundary layer mechanics. 

Fromi the concept of compressibility, we are just now 
evolving a comprehensive mechanics of supersonic flow 
to cope with the engineering needs of modern flight. 


Presented at the Newton Tercentenary Celebration, sponsored 
by the Royal Society, Cambridge University, England, July, 
1946. 

* Chairman. 
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Newton's ideas are as old as reason and as new as re- — 
search. Their timelessness is shown by three simple in- 
stances. 

Newton clearly stated that his laws of motion imply 
that it does not matter whether a body be fixed and the 
fluid flow over it or whether the body move through a 
fluid at rest. This kind of relativity is the justification 
for the wind-tunnel testing by the Wright Brothers and 
their followers in aerodynamic research. 

A second example of old and new is the modern 
wing theory which accounts for lift as the reaction to 
the force that communicates downward momentum 
continuously to the air. Froude applied the same 
reasoning to marine propulsion. 

The third instance is the newest style of all—jet pro- 
pulsion. Jet propulsion as a device of mechanics is cer- 
tainly Newtonian in principle. The propulsion is the 
third law reaction to the second law force, as measured 
by the rate at which momentum is created in the jet. 
While Newton is said to have speculated about a jet- 
propelled steam carriage, certain marine creatures were 
making use of a propulsive jet very much earlier and 
might claim priority in geological time. 


Two BRANCHES OF FLUID MECHANICS 


Leaving aside compressible and supersonic flows, I 
would like next to outline how the two great branches of 
fluid mechanics, frictionless flow and viscous flow, each 
starting with Newton, have independently developed 
over three centuries and now, in our time, have at last 
been brought together as a powerful instrument of 
aeronautical science. And finally, I shall have the honor 
to describe the successful engineering application of this 
new instrument of ancient lineage in the laboratory of 
the National Advisory Committee for Aeronautics; a 
laboratory named for Samuel P. Langley, an astrono- 
mer turned aeronautical scientist. 

I refer to the so-called laminar flow or low-drag wing, 
derived from theory and applied during the late war to 
both fighter and bomber airplanes. Its practical form 
is no longer a secret as, unfortunately, our enemies ob- 
tained samples from Mustangs that fell to their guns. 
This wing permits an airplane to fly faster and farther 
by virtue of its lower drag. By a coincidence, the pres- 
sure distribution that promotes smooth airflow also per- 
mits a higher flight speed before the onset of compressi- 
bility shock. 

To trace the pedigree of the low-drag wing, I must 
start with Newton, as I have intimated. 
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Frictionless Flow 


Newton’s dynamical concepts and his tactic of postu- 
lating an ideal fluid were taken up promptly by Daniel 
Bernoulli, who deduced the relation between pressure 
and velocity in his famous law of conservation of vis 
viva. d’Alembert further developed the concept of the 
frictionless fluid, showed how fluid should close in 
smoothly behind a body moving through it and so 
arrived at his Paradox of Zero Resistance. Euler 
formulated the mathematical statement of the problem 
and introduced a velocity potential such that its deriva- 
tive in any direction was the velocity component. 

These three were the founders of classical hydrody- 
namics. They were followed by the 19th century giants, 
Helmholtz, Kirchoff, Kelvin, and Rayleigh who con- 
tributed an understanding of the vortex, surfaces of 
discontinuity, and jets. Rankine showed how com- 
binations of sources and sinks could be used to give ideal 
forms for ships. 

At the beginning of this century, an elaborate struc- 
ture of applied mathematics had been built up from the 
hydromechanics of the ideal fluid but its conclusions 
were almost as remote from the behavior of real fluids 
as in Newton’s time. The subject was practically 
sterile. 

Such a static situation, however, was destined to be 
saved by the pioneers of early flight. It was discovered 
by practical men with no knowledge of hydromechanics 
that an arched wing profile was advantageous. Lan- 
chester in England, Kutta in Germany, and Joukowski 
in Russia independently accounted for this fact by the 
concept of circulation ardund the wing, increasing 
velocity over the upper surface. Even in the ideal fluid, 
which gave no resistance, a strong normal force or lift 
was predicted. Lanchester furthermore showed how a 
wing would cause the air to rise to meet it and then de- 
scend and that a vortex would trail behind each wing 
tip. He gave a physical picture of the flow which 
accounted for both lift and drag in a frictionless fluid 
and thus finally disposed of d’Alembert’s paradox of no 
drag. Joukowski, by a conformal transformation of a 
wing into a cylinder, computed the pressures on the wing 
from the velocity field. 

The picture was completed by Prandtl during the 
first World War. He paralleled and extended Lan- 
chester’s qualitative work and produced a complete 
airfoil theory of immense value to engineers. The 
theory was still in the domain of the ideal frictionless 
fluid and did not account for profile drag. This was left 
for wind-tunnel tests to determine. 


Viscous Flow 

Now to return to the other Newtonian concept, the 
viscous fluid. 

No such elaborate structure as hydromechanics was 
developed for the mechanics of viscous fluids. Experi- 
menters, or shall I say physicists, rather than mathema- 
ticians, dealt with such real fluids. However, early in 


the 19th century Stokes solved the equations of motion 
of a Newtonian viscous fluid and enunciated ideas that 
led to great advances. 

Osborne Reynolds, from careful observation of the 
flow of viscous liquid in tubes, distinguished two forms 
of motion, laminar and turbulent, with an abrupt transi- 
tion from the first to the second. By analysis, he found 
the determining factor to be a relation between velocity, 
tube diameter, density, and viscosity. This relation isa 
nondimensional coefficient, now known as the Reynolds 
Number, and is a criterion of dynamical similitude. It 
expresses the ratio between inertial and viscous forces in 
the fluid. 

That turbulent flow in pipes involved much more 
friction than laminar flow was recognized almost as soon 
as the controlling variables were distinguished. How- 
ever, the flow was much too complex to be amenable to 
simple analysis as for laminar flow, and experiments 
were confused by the indefinite nature of surface rough- 
ness. 

Prandtl in 1904 had shown that, at sufficient speed, 
the effects of viscosity must be largely confined to re- 
gions of fluid close to the boundary surfaces. This 
boundary layer contains the mechanism for transmis- 
sion of shear stress and profile drag to the wing. 

Blasius (1908) calculated, and verified experi- 
mentally, the force-momentum relations for the laminar 
boundary layer in terms of the velocity gradient or rate 
of shear in the Newtonian manner. Later, Tollmein 
(1935) discussed the stability of such laminar flow. 

Consideration of this theoretical work, supported by 
experimental evidence from B. M. Jones (1938) and 
others, led to the conclusion that a smooth surface and a 
falling pressure gradient are necessary to maintain an 
extensive laminar flow at high Reynolds Number. 
Otherwise, there is a prompt transition to turbulence. 

The force-momentum relations for a turbulent boun- 
dary layer on a flat plate, again following Newtonian 
dynamics, were derived independently by von Kar- 
man (1921) and Prandtl (1927). Subsequently, studies 
by G. I. Taylor, Dryden, and their collaborators 
accounted for skin friction and boundary-layer separa- 
tion with reasonable success. 

The outstanding practical characteristic of laminar 
flow was immediately evident—this is, that the drag is 
only a fraction of that for turbulent flow (Fig. 1). 
There was great engineering interest in means to prevent 
the ever present tendency toward instability and abrupt 
transition to turbulent flow. 

Theory showed that a falling pressure gradient was 
required to preserve laminar flow, and airfoils: of ar- 
bitrary shape were designed to effect such a pressure 
distribution. The results were disappointing. Wind- 
tunnel tests showed transition to turbulent flow at a 
Reynolds Number far too low for extensive laminar 
flow on full-scale airplane wings. Furthermore, transi- 
tion was sensitive to surface roughness, to any initial 
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Fic. 1. Skin-friction drag of flat plates. 


turbulence of the air stream, and even to acoustic vibra- 
tion. Although hope of success was slight, substantial 
drag reductiou could be obtained in no other way, since 
the profile drag of good wings at low incidence was at- 
tributable almost entirely to skin friction. 

Not only were wind-tunnel tests disappointing in 
1929-1930 but could be misleading as well. It was 
realized that all wind tunnels provided air streams of 
varying degrees of turbulence. An international ex- 
change of models took place, and much progress was 
made in explaining the discrepancy between tests of the 
same model in different wind tunnels on the basis of 
their inherent turbulence. 

This work probably stimulated the outstanding in- 
vestigations of G. I. Taylor and B. M. Jones in England 
and Hugh Dryden and E. N. Jacobs in America on wind- 
tunnel turbulence and its control. To Dryden, in 
particular, we are indebted for hot-wire anemometer 
equipment of beautiful precision and for research that 
explained fully the effect of screens in reducing wind- 
tunnel turbulence. 


Low-Turbulence Pressure Wind Tunnel 


In 1938 a low-turbulence pressure wind tunnel was 
built at the Langley Laboratory to test two-dimensional 
airfoils at a Reynolds Number corresponding to full 
scale. This tunnel had an area contraction of 20 to 1 
and a system of screens with the mesh and wire so 
selected that the wake from the screen did not contain 
vortices. Consequently, the screen itself did not intro- 
duce fine turbulence. In June, 1938, the first test of an 
airfoil designed to promote laminar flow indicated a drag 
coefficient of 0.0033 or about half of the lowest ever 
before recorded for an airfoil of comparable thickness. 
The turbulence measured was only a few hundredths 
of 1 per cent of the wind speed. 

Jacobs, with Dryden’s collaboration, was responsible 
for this wind tunnel and for what I believe to be a mile- 
stone in the progress of aerodynamics. The Langley 
research team now had a tool with which to untangle 
the confused state of our knowledge of viscous flow. 
They were greatly encouraged by B. M. Jones’ Wright 
Brothers Lecture, also in 1938, in which he clearly 
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showed from his experiments in flight that, by taking 
pains, laminar flow could be obtained for a consider- 
able distance behind the leading edge of a smooth 
wing. 

By May, 1941, the technique of airfoil design and 
testing at Langley had advanced to the point where low 
drags could be obtained for Reynolds Numbers corre- 
sponding to the flight of a medium bomber. However, 
the road was not easy and full success in the engi- | 
neering sense could-not have been attained by empiri- 
cal methods alone. Guidance of a general theory was 
required. 


DEVELOPMENT OF CLASSICAL HyDRODYNAMICS 


I now leave the story of viscous flow and return to the 
parallel development of classical hydrodynamics, which 
also stems from Newton. Recall that modern airfoil 
theory was well established by 1930 on the basis of an 
ideal fluid with circulation. Mathematical methods 
were available to compute the pressure distribution on 
certain types of airfoils which, however, had unsatis- 
factory characteristics and were not used. For all prac- 
tical purposes, airfoil theory in 1930 and practical air- 
foil design were unmarried. 

Analysis of pressure-distribution measurements led 
Jacobs to separate the effects of camber and angle of 
attack on lift and to assume them to be independent. 
Theodoresen showed (N.A.C.A. Report No. 383, “On 
the Theory of Wing Sections with Particular Reference 
to the Lift Distribution’’) that such superposition was 
justified by the theory of thin airfoils as stated by 
Munk and Glauert. As a result of this theoretical - 
support, a new family of N.A.C.A. airfoils was derived 
by Jacobs and Pinkerton in 1930. An arbitrary thick- 
ness distribution was used to form a series of sym- 
metrical sections. These were given various types of 
camber to form families of cambered airfoils. The 
series proved highly successful and was widely used 
throughout the world. A number of variations of 
thickness distribution gave further airfoil forms. This 
work, based on airfoil theory but essentially empirical, 
provided the basis for most of the airfoil work, especi- 
ally in Germany, Italy, and Japan up to and during the 
War. . 

The futility of efforts to obtain significant improve- 
ments by empirical variation of thickness distribution 
became apparent by 1935. What was wanted was 
stable laminar flow to a greater extent. To obtain 
this, complete control of pressure gradients was 
needed. 

In 1931 Theodorsen (N.A.C.A. Report No. 411, 
“Theory of Wing Sections of Arbitrary shape’) had 
published a theory of thick airfoils by means of which 
the pressure distribution for an arbitrary shape could be 
computed. The method used an inverse Joukowski 
transformation to transform an arbitrary airfoil to a 
nearly circular closed curve and, thence, to the corre- 
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ponding circle where the flow and pressures could be 
computed. Points on this curve are related to corres- 
ponding points on the circle by small angular and radial 
displacements. The pressure distribution as ordinarily 
computed from this theory was found to be not suffi- 
ciently precise near the leading edge for the prediction 
of local pressure gradients. Nevertheless, the theory 
was useful in showing how to obtain the favorable pres- 
sure drop over the forward portion of an airfoil and the 
sharp pressure recovery over the rear portion needed to 
prevent separation. Using approximate methods, air- 
foils of the N.A.C.A. families 2 to 5 were derived and 
tested. These families showed real advantages but are 
now considered obsolete. 

Further experience with the Theodorsen transforma- 
tion led Jacobs to the discovery that the airfoil, de- 
signed to give a desired pressure distribution, could be 
obtained with precision by an angular distortion of the 
curve into which the airfoil was transformed and then 
by correctiong both the airfoil and the corresponding 
pressure distribution by successive approximations. 
The method was systematized by Theodorsen and von 
Doenhoff in 1944 (N.A.C.A. No. L4G05, ‘‘Airfoil-Con- 
tour Modifications Based on ¢ Curve Method of Calcu- 
lating Pressure Distribution’’) and reduced to engineer- 
ing practice by Abbot, von Doenhoff and Stivers in 1945 
in their comprehensive ‘““Summary of Airfoil Daia”’ 
(N.A.C.A. No. L5C05). 

Engineering Applications 

Now a word as to what has been done with this theory 
and the low-turbulence wind tunnel in the way of appli- 
cation. This brings us into the realm of engineering and 
the art of design for a useful purpose. 
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Fic. 2. Pressure distribution. 


Fig. 2 shows the pressure distribution for a typical 
low-drag airfoil of the NACA 6 series as computed and 
as measured in the low-turbulence wind tunnel. Note 
the steadily falling pressure for the forward 60 per cent 
of the chord and the sharp pressure recovery after the 
position of minimum pressure. 


The greatest possible extent of laminar flow is to the 
position of minimum pressure. Beyond this point, tur- 
bulence is sure to occur. Transition from laminar to 
turbulent flow can, however, take place forward of the 
position of minimum pressure, depending on Reynolds 
Number, roughness, and possibly other factors. More 
research is required to clarify such questions. 

The position of minimum pressure is not far from the 
position of maximum thickness. This position deter- 
mines minimum drag by determining the extent of 
laminar flow. , 


In Fig. 3 are shown the lift and drag coefficients for a 


good conventional airfoil and for two low-drag sections 
of the NACA 6 series, all of the same thickness. The 
curves are faired through points representing experi- 
mental determinations in the low-turbulence wind 
tunnel. The reduction in drag effected is of the order of 
40 per cent. A method of computing such curves of lift 
and drag is presented in N.A.C.A. Report No. 4B05 
(‘A Concise Theoretical Method for Profile-Drag ,Cal- 
culation’’) and No. 4112 (‘‘The Theoretical Calculation 
of Airfoil Section Coefficients at Large Reynolds 
Numbers’) together with a comparison with wind- 
tunnel tests. Excellent agreement is shown. 
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Fic. 3. Experimental drag polars. 


The conventional airfoil shows a smooth variation of 
drag with increasing lift or angle of incidence. This is 
because the minimum pressure point moves slowly as 
the lift or angle changes. The sharp rise in drag for the 
NACA 6 series airfoils at two critical values of the lift is 
is caused by a sudden shift forward of the minimum 
pressure point and a consequent more forward transi- 
tion. 

This shift in the point of transition causes the charac- 
teristic ‘‘bucket’’ in the drag curve. It is, of course, the 
designer’s wish to widen this bucket to ensure low drag 
for a useful range of flight conditions. 

Fig. 4 illustrates the engineer’s close control of the 
characteristics of design. Here are the pressure dis- 
tributions for a series of airfoils designed to hold laminar 
flow for 30, 40, 50, or 60 per cent of the chord. The en- 
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Fic. 4. Theoretical pressure distributions. 

gineer may choose the one that best suits his operating 
conditions. He will recall that laminar flow is more 
difficult to maintain in proportion as the point of mini- 
mum pressure is farther back. He should be conserva- 
tive if he has reason to fear his airplane wing may not re- 
main smooth or may be required to operate under icing 
conditions in service. 

The present position is that, for subsonic flight, the 
airplane designer now has a remarkably powerful theory 
to guide his judgment. However, the theory has been 
developed without regard to compressibility and there- 
fore, does not apply at extremely high speeds. By a 
fortunate circumstance, the minimum pressure peak for 
these low-drag airfoils is lower than that for conven- 
tional airfoils, and, consequently, the critical flight speed 
at which sonic velocity is reached locally is somewhat 
higher. The great advantage of a combination of low 
drag with high critical speed was realized by the early 
NACA 16 series sections which are still used for pro- 
pellers where good characteristics at high speed are im- 
perative. 


CONCLUSIONS 


Let us pause to consider where we stand. Classical 
hydromechanics, still holding to its hypothesis of a 
frictionless fluid, predicts the normal pressures on the 
airfoil. The mechanics of a viscous fluid predicts the 
tangential forces on the airfoil. We assume, and the 
assumption is justified by experiment, that pressures 
due to the external flow are transmitted unchanged 
through the boundary layer. Thus, at last, we have the 
marriage of the two main branches of fluid mechanics 
which started with Newton’s speculations on the nature 
of fluid resistance in his two imaginary fluids. 

I should report for the record that the N.A.C.A. re- 
search that culminated after some 15 years of effort in 
this important advance was planned, guided, and in- 
spired by Dr. George W. Lewis, Director of Aero- 
nautical Research. His research teams included Jacobs, 
Theodorsen, Dryden, Abbott, Stack, and many others 
whose names will be found on N.A.C.A. reports since 
1930. 

While emphasizing this recent advance in fluid 
mechanics, I should be less than frank if I neglect to 
mention that there is a great deal more to learn. Speci- 
fically, the prediction of transition and separation of the 
boundary layer is not yet possible, and the entire 
transonic régime is a mystery. When the speed of 
flight equals the speed of sound, the flow is part sub- 
sonic and part supersonic. There is no adequate theory 
for such mixed flows. Furthermore, the wind tunnel, 
on which engineers have always depended for design 
data when theory would not serve, completely fails 
them at the speed of sound. Shock waves from the 
model are reflected back by the tunnel walls. 

This surely is a challenging opportunity for research. 
Let us hope it will not take 300 years. 





Letters to 


Dear Sir: 

The method of Ramberg and Osgood! of describing stress- 
strain curves by three parameters provides a convenient method 
of specifying the proportional limit of a material. , 

The importance of the proportional limit in compressive in- 
stability calculations has always been realized, but, because of 
the difficulty of precisely determining this property experiment- 
ally, it is usually neglected in specifying material properties. 

The three-parameter method obviates experimental deter- 
mination of the proportional limit, since the stress-strain curve 
is completely defined by an analytical expression, the constants 
of which are experimentally ascertained. If we accept the defi- 
nition of proportional limit as the stress at which the strain has 
departed 0.01 per cent from that obtained from elasticity, it is 
a simple matter to derive an analytical expression for the pro- 
portional limit. 

In the body of the paper,! the following expression for a stress- 
Strain curve is given when m, = 0.7. 


Ee/S; = o + 3/70” (1) 


the Editor 
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From Fig. 2, it can be seen that 
én = 0.0001 + (S,:/E£) (2) 
Combining Eqs. (1) and (2) and solving for the proportional limit 
stress ratio, 
on = Sp/S; = [2.833 X 10-4 (E/S,)]/” (3) 


This expression is depicted on Fig. 1, on which curves are drawn 
in terms of the three basic parameters. 


REFERENCE 


1 Ramberg, W., and Osgood, W. R., Description of Stress-Strain Curves by 
Three Parameters, N.A.C.A. T.N. No. 902, July, 1943. Osgood, W. R., 
Stress-Strain Formulas, Journal of the Aeronautical Sciences, Vol. 13, No. 1, 
p. 45, January, 1946, 


GEORGE GERARD 
Principal Research Engineer 
Republic Aviation Corporation 





Dear Sir: 

To facilitate flutter computations in the design stages of an 
airplane and to enable the flutter analyst to gain a better insight 
into the critical and noncritical parameters of the design it is 








TABLE 2 
i 0.25 0.50 1.00 1.25 1.67 2.00 2.50 
+2.0 +2.5 +2.0 +1.5 +0.75 -—0.5 +0.1 
% error) 
( —-3.9 —4.4 —2.7 —1.4 +1.4 +3.5 +7.5 


— 
(% error) 





suggested that finite, yet close, approximating functions be sub- 
stituted for the commonly used tabular values of Theodorsen’s F 
and G aerodynamic functions. 

For a range of V =v/bw from V = 0 to V = o, it is suggested 
that the following approximating functions of V be used: 


F @ 1/.{1 + tanh (V/10)} 
Ga —V/8e— V/s 


It is to be observed that these approximatiors yield values 
rigorously in accord with the true values of F and G for the 
limiting conditions V = Oand V = o. 

For intermediate values of V, the closeness of approximation 
may be ascertained by a study of Table 1, 

Despite the closeness of fit of the approximating functions, as 
seen in Table 1, they are, nevertheless, still too complicated to 
allow the flutter determinant to be solved explicitly for V. 
Therefore, the following further approximations are suggested: 


Fa '/,(1 + V/10) 
Ga —V/8(1 — V/4 + V?/32) 


For the range of significant V’s encountered in the analysis of 
control surfaces of transport airplanes (V = 0.25 to V = 2.50), 
these latter approximations hold as closely as the former did over 
a wider range (V = 0.25 to V = 5.00). Closeness of fit may be 
judged from Table 2. 

Even these errors are not so serious as one would at first sup- 
pose because of the manner in which F and G enter the total 
aerodynamic vector. For example, consider Theodorsen’s 
App = Rog + tlpg. For it, the real and imaginary parts may be 
written as 


Rog = —Bgi — Q(c)(2G)V + [Bes + Q’(c)(2F)]V? 
Ib = [Bg2 + Q(c)(2F)]V + Q’(c)(2G) V? 


where all the symbols are Theodorsen’s save for V and the func- 
tions of c, Q(c) and Q’(c). The Q functions are derived in the 
process of simplification of Rpg and Jpg. All of the aerodynamic 
vectors may be written in exactly similar form. 

For even rougher approximation and correspondingly greater 
simplification of the aerodynamie terms, (2G) may be taken as 
zero and (2F) as unity. The real parts then become quadratic 
and the imaginary parts become linear. This corresponds ex- 
actly with experience in plotting the aerodynamic terms using the 
true values of F and G and when the aerodynamic terms were 
plotted for a workable range of V = 0 to V = 2.50. 

It is the hope of the writer to elaborate on the last paragraph 
above in an article at some future date. 

W. RAYMOND GRIFFIN 
Research Engineer 














TABLE 1 
V 0.25 0.50 1.00 1.25 1.67 2.00 2.50 3.33 5.00 
0.5037 0.5129 0.5395 0.5541 0.5788 0.6030 0.6245 0.6650 0.7276 
(true) 
0.5120 0.5250 0.5500 0.5625 0.5835 0.6000 0.6250 0.6592 0.7310 
(appr ox.) 
Error +1.5 +2.5 +2.0 +1.5 +1.0 —0.5 0 —1.0 +0.5 
(%) (approx.) 
(-—G) 0.0305 0.0577 0.1003 0.1165 0.1378 0.1510 0.1660 0.1800 0.1886 
(true) 
(-—G) 0.0293 0.0552 0.0974 0.1144 0.1374 0.1516 0.1672 0.1810 0 1790 
(approx.) , 
Error —4.0 —4.4 —3.0 —-1.8 —0.3 +0.5 +0.9 +0.5 —5.1 


(%) 
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Superaerodynamics, Mechanics of Rarefied 
Gases 


HSUE-SHEN TSIEN* 
California Institute of Technology 


INTRODUCTION 


Se IN 1934 published an article on the aerody- 
namics of highly rarefied gases, a branch of fluid 
mechanics which he called superaerodynamics. At that 
time, however, with the means of propulsion then avail- 
able, flight at extreme altitudes did not seem to be 
realizable. Therefore, superaerodynamics has been 
considered as a subject of academic interest rather than 
one of practical engineering importance. With the re- 
cent perfection of the rocket as a propulsive power plant, 
the situation is radically changed and there should be 
no limit to the altitude that can be reached by an air- 
craft. There are even indications that the optimum 
flight altitude of long-range rocket airplanes is approxi- 
mately 60 miles. At these and higher altitudes, the air 
density is so low that the fluid must be thought of as one 
having a coarse structure and not as a continuous 
medium of conventional fluid mechanics. Then con- 
cepts of superaerodynamics are needed to guide the 
design of such an aircraft. 

Besides its applications to high-altitude flights, 
superaerodynamics should have many applications to 
industrial processes in which low density gases are in- 
volved. The knowledge of superaerodynamics should 
be of invaluable help in improving the efficiency of 
these processes. For instance, the pumping of low 
density gases will be of more and more importance 
because of the increasingly frequent use of high vacuum 
in distillation and other chemical engineering processes. 
The improvement in the design of such gas pumps defi- 
nitely needs the understanding of the principles of 
mechanics of low density gases. It is the purpose of 
the present paper to discuss the fundamental concepts 
of this new branch of fluid mechanics and to indicate 
some of the results already obtained. The field of 
superaerodynamics is still relatively uncultivated. 
Here and there the ground has been broken by a few 
physicists. However, further effort is required to de- 
velop this branch of fluid mechanics into an aid in 
engineering design and research. 


MEAN FREE PATH AND REALMS OF FLUID MECHANICS 


As a first approximation, the gas may be considered 
as an aggregate of rapidly moving particles which are 
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constantly colliding with each other. The influence 
of the particles on each other can be conveniently neg- 
lected until they are so close together that a ‘‘collision’ 
takes place. Then the coarseness of the structure of the 
gaseous medium can be expressed by the parameter / 
which is the distance the particles or molecules travel 
between collisions. Since the instantaneous velocity 
distribution and density distribution in the gas are far 
from uniform, one can only conveniently use the statis- 
tical average of a quantity instead of the instantaneous 
value of the quantity. The distance / is then the 
statistical average over the billions and billions of mole- 
cules concerned. This average / is called the mean 
free path of the gas. If the mean free path is small in 
comparison with the dimension of the flow field or the 
dimension of the body in the flow field, then the gas 
can be considered as a continuum and the ordinary gas 
dynamics is sufficient for the, analysis of these flows. 
If the mean free path is not negligible when compared 
with the dimension of the body, then the effects of the 
discrete character of the gas must be taken into ac- 
count in the calculation. Then if L is the linear dimen- 
sion of the body, superaerodynamics can be defined as 
the aerodynamics of flows where the ratio //L is not 
negligible. 

The free mean path / can be calculated from the mean 
velocity of the molecules, the density of the gas, and the 
“effective radius of the molecules.”’ However, there is 
no direct way of measuring the effective radius of the 
molecules; therefore, it is more fruitful to express the 
quantity / in terms of a measurable quantity such as 
the viscosity of the gas. This is easily done by the 
following consideration. If the gas flows in the x- 
direction with the macroscopic velocity u(y) (Fig. 1), 
the gradient of velocity is then 0u/Oy. Now the gas 
molecules in a lower layer move into an upper layer with 
the average velocity 7 of the molecules. The distance 
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Fic. 1. Shearing action in parallel flows. 
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they will travel before they lose their identity by col- 
lision with other molecules is/. The difference in macro- 
scopic velocity of the layer where the molecules origi- 
nated and the layer where the molecules are mixed is 
then 1(0u/Oy). The mass of molecules crossing a 
unit area of the layer proportional to pv where p is the 
density of the fluid. Therefore, momentum exchange 
between layers is proportional to pi/(Ou/dy). This is 
the viscous shearing stress. Since the coefficient of 
viscosity is defined as 


r = p(du/dy) 


» is proportional to pvl. The proportionality constant 
can be calculated by the kinetic theory. The most 
accurate calculation is that due to Chapman? who gives 


Bh = 0.499pol (1) 


The average velocity 0 is closely connected with the 
“velocity of sound”’ a by the equation 


8 |p 8 
v= YE V = —_ a (2) 


where y is the ratio of specific heats. From Eqs. (1) 
and (2), the mean free path is given in terms of the 
kinematic viscosity v = u/p and the vefocity of sound a 
by 


L = 1.255 V y(v/a) (3) 


For air, the quantity /(p/po) in inches is tabulated in 
Table 1, where p and fp are the pressure and standard 
pressure (one atmosphere), respectively. Thus, for 
ordinary pressures, the mean free path is truly negligible 
in comparison with the body dimensions, and therefore 
gas dynamics is sufficient for aerodynamic calculations. 
However, at extremely high altitudes where the pres- 
sure may be only one-millionth of the pressure at the 
surface of the earth, the mean free path will be com- 
parable to the body dimensions. A clear demonstration 
of this fact is perhaps that of Fig. 2 where the mean free 
path calculated by Maris? is plotted against the altitude. 
Maris assumed that the air temperature above 80 km. 





TABLE 1 
Mean Free Path for Air 





Tem- Tem- 
perature, perature, 

< l(p/po), In. “‘. l(p/po), In. 
0 2.32 X 10-* 260 5.33 X 1076 
20 2.55 X 10-° 280 5.57 X 1078 
40 2.78 X 10-* 300 5.85 X 107-6 
60 3.00 X 107° 320 6.05 XK 10-6 
80 3.23 X 10-¢ 340 6.28 <x 107° 
100 3.46 X 10-6 360 6.52 XK 10-6 
120 3.69 X 10-6 380 6.77 X 107° 
140 3.92 X 107° 400 7.00 X 10-6 
160 4.16 X 107-6 420 7.24 X 107-6 
180 4.40 X 1076 440 7.49 X 107-6 
200 4.62 X 107-6 460 7.73 X 10-6 
220 4.85 x 10-6 480 7.96 X 107-8 
240 5.10 K 1076 500 8.20 K 10-6 





po = Standard pressure, one atmosphere. 
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Fic. 2. Mean free path at different altitudes; SD = summer 
day, SN = summer night, WD = winter day, and WN = winter 
night (H. B. Maris). 


(50 miles) is constant and has the value 97°C., —43°C., 
—23°C., and —53°C. for summer day, summer night, 


winter day, and winter night, respectively. At 50 miles_ 


altitude, the mean free path is seen to be nearly 1 in. 
The air there is certainly not a continuum. 

The flow conditions near the wall for the case in which 
the mean free path / is small but not negligible compared 
with the body dimension LZ or the thickness 6 of the 
boundary layer were first investigated theoretically by 
J. C. Maxwell in 1879. It was found that the gas no 
longer sticks to the surface but slips over the surface 
with a definite velocity. This phenomenon was also 
suggested by the experiments of Kundt and Warburg 
in 1875, on the damping of a vibrating disc by gas at 
low pressure. Later experiments and theoretical study 
completely verified this conclusion. This type of flow 
can be called the slip flow. Because of Eq. (3), 


1 LM 


Vv 
mt tot (4) 


l 

_* 
where M is the Mach Number of the free stream and 
Re is the Reynolds Number of the flow referred to the 
length L of the body. For extremely small Reynolds 
Number, L/6 ~ 1, then 


l/i~ M/Re, Re<1 (5) 


For large Reynolds Number, it is well known that 
L/i ~ V Re. Then, 


1/5 ~ M/V Re, Re>1 (6) 


Therefore, if the interval 1/100 < 1/5 < 1isconsidered 
as the proper range for the slip flow, then in the plane 
of M and Re, this realm of fluid mechanics occupies a 
region as shown in Fig. 3. The region below this slip- 
flow region belongs to the realm of the conventional 
gas dynamics with the usual boundary conditions at the 
wall. 

If the mean free path is much larger compared with 
the body dimensions, one enters an entirely new realm 





of fl 
mol 
chat 
the 


‘of fc 


of n 
dete 
streé 
of tl 
the 
plific 
sidet 
to t 
mole 
can 

If 
free 
occu 


This 
mole 
mec] 
the | 
impc 
and 

offer 
itself 
still 

Rea 


STR 


Th 
are tl 
the e 
veloc 
press 
tivel) 


In th 
stress 
vecto 





"i 





LAJ 





su Ym Vee ewes hOlUlUrlU 





SUPERAERODYNAMICS 655 


of fluid mechanics. Here the chances for the collision of 
molecules among themselves are much smaller than the 
chances for the collision of molecules with the wall or 
the surface of the body. Therefore, for the calculation 
of forces, one need only consider the impact of a stream 
of molecules with a velocity and energy distribution 
determined by the thermal equilibrium in the free 
stream, the Maxwellian distribution. The re-emission 
of the molecules from the surface will be governed by 
the accommodation coefficient. But the greatest sim- 
plification comes from the fact that one need not con- 
sider the distortion of the Maxwellian distribution due 
to the collision of the re-emitted molecules with the 
molecules in the stream. The realm of fluid mechanics 
can thus be called the free molecule flow. 

If one takes the limiting ratio of 1//L to be 10 for the 
free molecule flow, then this realm of fluid mechanics 
occupies a region of the M-Re plan given by 


This is shown in Fig. 3. The region between the free 
molecule flow and the slip flow is the realm of fluid 
mechanics where the collision between molecules and 
the collision of molecules with the wall are of equal 
importance. The problem is extremely complicated 
and no satisfactory theoretical solution can yet be 
offered. However, it is certain that as far as the fluid 
itself is concerned, the characteristic parameters are 
still the Mach Number M and the Reynolds Number 
Re as used in gas dynamics. 


STRESSES AND BOUNDARY CONDITIONS IN SLIP FLow 


The fundamental equations governing the fluid flows 
are the continuity equation, the dynamic equations, and 
the energy equation. If m4, “2, u3 are the components of 
velocity in the x), x2, x3 directions, and p, p, T the 
pressure, density, and temperature of the gas, respec- 
tively, these equations are* 


(Op/Ot) + (Opu;/Ox;) = 0 (7) 
Ou; Ou; Op Ori; > 
— Sete ae yes nm 8) 
Oe tems Ox; Ox, AWM, { 
( o oO A: ot i op Og: _ 
tt = “adie ) ~ Of ” Oxy 


o” 
ag (Tijt;) (9) 


In these equations, the r’s are the components of the 
stress tensor, the g’s are the components of the heat flux 
vector, and c, is the specific heat at constant pressure. 


* The summation convention is used here. For instance, 


Ou K Ou; Ol2 Ou; 


OX, Ox) OX» OX; 
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3 77 
U0 Re 
Fic. 3. Realms of fluid mechanics. 


For gases under ordinary pressures, the stresses are 
the products of the coefficient of viscosity and linear 
combinations of the first order derivatives of the veloc- 
ity components. The components of heat flux are 
simply the products of coefficient of heat conduction 
and the components of the gradient of tempera- 
ture. 

These relations between the stresses and the heat flux 
with the flow quantities are confirmed by the kinetic 
theory of nonuniform gases developed by Maxwell, 
Boltzmann, Chapman, D. Enskog, and others. In this 
theory, the distribution function of the molecules in the 
velocity space on the probability of molecules having a 
particular range of velocities plays the predominating 
role. 

If, in the macroscopic sense, the gas is at rest 
with uniform temperature, the probability function is 
the well-known Maxwellian distribution. In a non- 
uniform gas, the first order correction to this Maxwel- 
lian distribution gives the ordinary viscous stresses 
and the heat flux by conduction, as observed by experi- 
ments. 

However, the theory 
correction to the Maxwellian distribution which intro- 
duces highly complicated additional terms to the 
stresses and the heat flux. These additional terms, 
while negligible for gases under ordinary pressure, are 
important for rarefied gases. In fact, they form the 
basis for interpretation of many experimentally ob- 
served phenomena, such as the operation of a radiom- 


also gives a second order 


eter. 
The investigation of Burnett‘ gives the following 
forms of the heat flux and the stresses: 


oT pe ou; OT pe? 2 o) (7281) rs 


m= XZ +h Toe, de, + TLS MA. Oe, 
; 2 2 2aT 
, OM <4 (+. va fe Nig legh ly 2f) ;) 
Ox; OX; p Ox; p Ox, pT Ox; 


(10) 
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2 Ou 2 
Ty = —Qyey + K,- —* ey + Ke 2 | - 2(} 2e = 


p Ox» p Ox; p Ox; 
2 D2 
Oe Oe 264 Oe) + eT 
Ox; OX, Oxr pT Ox,Ox; 
uw? Op OT wu? OT OT pao 
iar tee, oe, ac — 11 
; ppT Ox, Ox; : pI? Ox, Ox; + Ke p Cixlxy (11) 


In these equations, \ is the coefficient of heat conduc- 
tion and yu is the coefficient of viscosity. The first 
terms of the equations are thus the ordinary heat flux 
and viscous stresses. 61, 02, 63, 0, 06; and Ki, Ke, Ks, 
Ky, Ks, Keg are pure numbers. Furthermore, 


, Ou; ous) Ouz 
a 1/f St 4. 9} — t/, 3 12 
Ci At Se, /s ed ij ( ) 


where 6;; = 1ifi = j and 6,, = Oif7 #7. Any tensor 
A,,; with a bar over it has the following meaning: 


Ay _ 1/(A y+ Ax) fad /3 Andi (13) 


The approximate values for 6’s and K’s are 


a= (2 - Tex) 
/ 4\2 p aT] 


= 84 = -3,0=3,4= 3 (5-7) (14) 


p aT. 

ma$(2-T%), 
8\2 p aT 

3T du 


K, = 2,K,; = 3, K, = 0, K, = — —, Ke = 8 (15) 
wp aT 


More accurate values for K’s are given by Burnett‘ for 
molecules of rigid elastic spheres. They differ little 
from those of Eq. (15). 

It is immediately clear by examining Eqs. (10) and 
(11) that the additional terms in the heat flux and the 
stress tensor are only of importance in rarefied gases. 
As an example, the magnitude of the ratio of the second 
terms to the first terms of these equations will be esti- 
mated presently. Let U be the free stream velocity 
and 7) the free stream density and temperature, c, 
the specific heat at constant pressure. For gases, the 
Prandtl number c,u/d is rarely unity. Then these 
ratios have the same order of magnitude given by 





Mr, sho tteo y= 
L Poly! o poUL Cpl o R 


where M is the free stream Mach Number and Re the 
Reynolds Number. According to Eq..(4), the order of 
magnitude of the additional heat flux or stresses to the 
ordinary heat flux and stresses is then M(//L). That 
is, the additional terms are only of importance if the 
product of Mach Number and the ratio of mean free path 
to body dimension is large. Therefore, for the gas 


dynamical flows with small mean free path, the ordi- 
nary heat conduction terms and viscous stress terms 
suffice. This is also true for the slip flows if the Mach 


Number is small. For slip flows at high Mach Number, 
the conventional Navier-Stokes equations are no longer 
a true description of the physical relations and the 
complicated forms of heat flux g; and the stresses 7, 
given by Eqs. (10) and (11) must be used. 

Since the additional terms in gq, and 7,; contain higher 
derivatives than the first order, they will raise the order 
of the system of partial differential equations given by 
Eqs. (7), (8), and (9). Hence, to solve these equations, 
one needs more boundary conditions than in the case 
of gas dynamics. In fact, the question of proper bound- 
ary conditions for slip flows of large Mach Number 
where such additional terms are required, has not yet 
been answered. This is then one of the fundamental 
problems in superaerodynamics. The theoretical study 
of the boundary conditions must, of course, be based on 
the kinetic point of view. Parallel to such theoretical 
investigation, it seems advisable to obtain also an ex- 
perimental check on the validity of Eqs. (10) and (11), 
One of the possible experiments would be to observe 
the flow of rarefied gas between two concentric cylinders, 
one rotating and the other held fixed. If the length of 
the cylinders is large in comparison with the gap be- 
tween the cylinders, the flow will be two-dimensional 
with the radial distance as the only parameter. The 
physical situation is thus the simplest, and amenable 
to detailed examination. Indeed, the same experiment 
has been performed by Millikan and his collaborators 
at low rotative speeds for the same type of investigation. 
The high rotative speeds necessary to obtain high flow 
velocities will certainly introduce difficulties in the 
technique but these do not seem to be insurmountable. 

The results of the kinetic theory discussed in the 
previous paragraphs are based upon the assumption 
of simple molecules with the only intrinsic energy due 
to translation. Generally, however, the gas molecules 
have three kinds of intrinsic energies: the translational 
energy, the rotational energy, and the vibrational 
energy. When the gas is in equilibrium, the distribution 
of the total energy among the three forms is uniquely 
determined by the properties of the gas molecules con- 
cerned. If this equilibrium is destroyed by a sudden 
change in the external conditions, such as an expansion, 
new equilibrium will be reached by numerous collisions 
of the molecules under these new conditions. However, 
it is known through ultrasonic dispersion, etc., that the 
process of reaching equilibrium in the vibrational 
energy isa slow one. In other words, a large number of 
collisions is necessary to excite the vibration degrees of 
freedom of the molecules properly. Since the average 
velocity of the molecules is 3, the average distance 
traveled by the molecules per unit time is 7 1. The 
number of collisions made by the molecules per unit 
time is then 0/]. Therefore, the time ¢, generally called 
the relaxation time, necessary to excite the vibrational 
energy is proportional to the mean free path at a given 
temperature. In other words, the relaxation time is 
inversely proportional to the pressure of the gas. The 
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; TABLE 2 
Relaxation Time ¢ for the Excitation of Vibrational Energy at 
One Atmosphere (Reference 6) 











— Pure CO,.——_—~ O.— ~ 
Temperature, Temperature, 
a t X 10° sec. "S. t X 108 sec. 
19.2 10.8 40.7 4.4 
47.1 10.6 49.7 2.9 
58.1 8.1 77.9 4.0 
88.7 7.9 100.2 3.6 
138.0 3.1 
140.6 3.7 
141.4 3.3 





measured value of ¢ for several typical gases is given in 
Table 2 (reference 6). Therefore, at extremely low 
pressures, the relaxation time may be of the order of 1 
sec. If this is coupled with high flow velocity, the gas 
must have passed the body before the vibrational de- 
grees of freedom are excited. Therefore, in low density 
gas flows, the vibrational energy of the gas can be 
considered as fixed at.the “free stream’’ value. This 
decrease in the effective number of degrees of freedom 
of the gas molecules tends to raise the value of y, the 
ratio of specific heats. It is fortunate that for air at 
ordinary temperatures, the equilibrium vibrational 
energy is small, therefore this freezing of vibrational 
degrees will not greatly alter the value of y. However, 
at higher temperatures and for other polyatomic gases, 
the effect may be appreciable. 


BOUNDARY CONDITIONS FOR SLIP FLOWS OF SMALL 
Macu NuMBER 


As concluded in the preceding section, the slip flows 
of small Mach Number can be calculated by using the 
Navier-Stokes equations. The proper boundary condi- 
tions for such flows were first investigated by Maxwell. 
Millikan’ has reduced Maxwell’s consideration to a 
simple form. A more detailed theoretical treatment of 
the problem by Boltzmann’s H-theorem was given by 
Epstein.’ The result of these investigations shows that 
if the Mach Number of the flow is small, then the flow 
velocity along the wall, the normal gradient 0u/Ody of 
u, and the temperature gradient 07/0x along the wall 
are connected by the relation*® 


st?yta & 


— (16) 
B oy 4 pT Ox 


at y = 0, where £ is the coefficient of sliding or external 
friction and yu the coefficient of viscosity. The second 
term in Eq. (16) is the velocity due to temperature gradi- 
ent and is called the creep velocity. The ratio of u and 
8 is given by 


u/B = 0.998 [(2/f) — 1]! (17) 


f is the fraction of the tangential momentum of the im- 
pinging molecules transmitted to the wall. Maxwell 
himself interpreted the fractional value of f as meaning 
that a fraction f of the surface reflects diffusely and the 











TABLE 3 
Values of f (Reference 5) 
f 
Air or CO, on machined brass or old shellac 1.00 
Air on mercury 1.00 
Air on oil 0.895 
CO, on oil 0.92 
Hydrogen on oil 0.925 
Air on glass 0.89 
Helium on oil 0.874 
Air on fresh shellac 0.79 





remainder specularly. The value of f is listed in Table 3 
as given by Millikan.’ It is seen from Eq. (17) that if 
l is small compared with the scale of the velocity 
gradient or the thickness 6 of the boundary layer, then 
8 is large and u at the wall is small. This is the reason 
why in gas dynamics one generally imposes the bound- 
ary condition « = 0 at the wall. If //é is not negligible, 
then there is a slip of the gas at the wall. 

Table 3 shows that the value of f is very nearly 
unity. Therefore, the molecules, after striking the 
surface, have a strong tendency to reflect diffusely. 
In the scale of molecular dimensions, the surface of a 
solid body must be extremely jagged even if it is highly 
polished in the macroscopic sense. Then one cannot 
expect the molecules to reflect in any uniform direction. 
The tendency toward diffusive reflection is thus, per- 
haps, anticipated. Furthermore, a closer study of the 
behavior of the molecules colliding with the solid sur- 
face seems to indicate the temporary adsorption of the 
molecules on the surface. This combined structure of 
gas molecule and the surface materials then breaks up 
with a certain average rate determined by the “‘life- 
time’’ of the combined structure. The molecules thus 
freed are re-emitted from the surface. With this pro- 
cedure, the molecules re-emerge from the surface with 
all their prior history obliterated. In particular, the 
direction of re-emission iscompletely dependent upon the 
direction of impinging. In other words, the molecules 
re-emit diffusely. 

Similar to this velocity discontinuity, there is a tem- 
perature discontinuity at the wall first found by Smolu- 
chowski in 1898. Namely, 


«(T — Ty) = MOT/dy), y =0 (18) 


where T is temperature of the gas, 7, is the tempera- 
ture of the wall, and x is the coefficient of temperature 
discontinuity. This coefficient x can be expressed in 
terms of the mean free path by a consideration similar 
to Maxwell’s.? The result is 


1 2-—a v 1 
-= .996( je —— 19 
K a /y+1nG (19) 





where a is the accommodation coefficient introduced by 
Knudsen.’® a@ can be defined as the fractional extent 
to which those molecules that fall on the surface and 
are reflected or re-emitted from it, have their mean 
energy adjusted or “accommodated” toward what it 





658 JOURNAL OF THE AERONAUTICAL SCIENCES-DECEMBER, 1946 


m Se i l 72 -_ 
—_ = 0.998 ( _ )é = 1.253 V ( — i) (22) 
BR f : 


would be if the returning molecules were issuing as a 
stream out of the mass of gas at the temperature of the 
wall. This adaptation of the impinging molecules to 
the wall conditions is, of course, expected from the 
temporary adsorption of molecules on the surface, as 
mentioned before. Thus, if £; is the energy brought up 
to unit area of the wall per second by the incident mole- 
cules, and E, the energy carried away by the re- 
emitted molecules, and £, the energy that would be 
carried away if the gas is at the wall temperature, then 


E, — E, = a(E; — Ey) (20) 


Of course, strictly speaking, a is different for the three 
different types of molecular energies: translational, 
rotational, and vibrational. However, experimental 
evidence seems to indicate that one coefficient is gen- 
erally sufficient. The older measured values of a 
are given in many books on kinetic theory of gases."! 
Wiedmann” recently determined the value of a for air 
on metals. The results are given in Table 4. Wiedmann 
concluded that the value of the accommodation coeffi- 
cient is independent of the nature of the technical sur- 
face of the metal. This again supports the concept of 
the microscopic character of the surface collision. 


TABLE 4 
Accommodation Coefficient a for Air (Reference 12) 
—-Accommodation Coefficient— 








Surface Description Minimum Maximum 
Flat black lacquer on bronze 0.881 0.894 
Bronze, polished 0.91 0.94 
Bronze, machined 0.89 0.93 
Bronze, etched 0.93 0.95 
Cast iron, polished 0.87 0.93 
Cast iron, machined 0.87 0.88 
Cast iron, etched 0.89 0.96 
Aluminum, polished 0.87 0.95 
Aluminum, machined 0.95 0.97 
Aluminum, etched (0.89) (0.97) 





S.trp Flows AT SMALL Macu NUMBERS 


If the Mach Number is extremely small, the effect 
of the compressibility of the gas cam be neglected and 
the fluid can be treated as incompressible. In order 
that the mean free path will be comparable with body 
dimension, the Reynolds Number must also be small. 
The flow can then be calculated by the well-known 
method of Stokes'* who neglected the inertia terms, 
which are small in comparison with the viscous and 
pressure terms. The resistance of a sphere of radius 
R moving with a velocity U through a fluid of viscosity 
and coefficient 6 of sliding friction is given by Basset.!4 
If the drag coefficient Cp is defined as the drag divided 
by (p/2) U*rR? then 


9 9 / 
Cp = r2 Be (21) 
Rell + 3(u/BR) 


where Re is the Reynolds Number UR/v. According to 
Eqs. (3) and (17), 


Qt 


R f Re 


Hence, Eq. (21) can be rewritten as 


12 1 + 2.506Vy [(2/f) — 1](M/Re)| 


\, 





> Re} + 3.759V 7 [(2/f) — 1](M/Re)\ 


(sphere) (23) 
Fig. 4 shows the drag coefficient as function of M and 
Re with y = 1.405 and f = 1. It is seen that there is 
an appreciable reduction in the drag coefficient due to 
the slip phenomenon. It also introduces a large effect 
of Mach Number, although this is not the compressi- 
bility effect in the ordinary sense. 
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Fic. 4. Drag coefficient Cp as functions of Mach Number M 
and Reynolds Number for spheres of radius R, f = 1, y = 1.405 


The similar case of the uniform two-dimensional mo- 
tion of a circular cylinder can be obtained by a slight 
modification of Lamb’s solution.'® If the drag coeffi- 
cient is defined as the drag per unit length divided by 
(p/2) U?(2R) where R is the radius of the cylinder, 
C> = 

4ar/ Re log = — 1.28107 + 1.253Vy (?- 1) zs 

Re “ont a f Re 
(cylinder) (24) 
where Re is again the Reynolds Number UR/y. The 
results of calculation are again plotted in terms of M 
and Re in Fig. 5. The reduction-of drag due to slip is 
again evident. 

For the fully established laminar flow through a 
circular tube of radius R, the velocity u as a function of 
the radial distance from the axis of the tube is deter 


mined by 
ldp _1 af A *) (25) 
p dx r dr dr 
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Fic. 5. Drag coefficient Cp as functions of Mach Number M and 


Reynolds Number for cylinders of radius R, f = 1, y = 1.405. 


The pressure gradient along the tube is a constant. 
Therefore, Eq. (25) yields 

du/dr = (1/2u)(dp/dx)r (26) 
and 

u = (1/4y)(dp/dx)r? + uo (27) 
where % is the velocity at the center of the tube. The 
boundary condition that 

Bu;-r = — u(du/dr),or 
then gives 

— dp/dx = 


(4u/R*)/[1 + 2(u4/BR)]u (28) 


The relation between the average velocity U and the 
maximum velocity uw can be easily obtained by integrat- 
ing Eq. (27). Thus 


Vill + 2(u/BR)]/['/2 + 2(u/B8R)]} (29) 
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Fic. 6. Pressure drop coefficient A as functions of Mach 
Number M and Reynolds Number Re for flow in a tube of radius 
R,f = 1,7 = 1.405. 
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As the ratio of mean free path to the radius increases, 
the value of u/8R will increase according to Eq. (22). 
Then Eq. (29) shows that the average velocity U ap- 
proaches the maximum velocity m. In other words, 
the velocity distribution across the tube becomes more 
uniform. This is, of course, expected from the “‘slip’”’ 
at the wall. 
If the pressure drop coefficient A is defined as 


A = (4R/pU*)(— dp/dx) (30) 
then Eqs. (19) and (20) give immediately 
64 1 


A== — 
Re 1 + 10.024 Vy [(2/f) — 1](M/Re) 





(31) 


where Re is the Reynolds Number based upon the diam- 
eter, U(2R)/v. Here again the effect of slip is to reduce 
the resistance coefficient. For extremely small mean 
free path the value of the second fraction is unity and 
the equation reduces to the well-known formula for 
Poiseuille flow. Fig. 6 is a graph for A as function of M 
and Re, assuming as before f = 1 and y = 1.405. Ex- 
periments on the slip flow in a tube were made by 
Knudsen,!*: 17 Gaede,!® and Ebert. 2° The correct- 
ness of Eq. (31) is fully verified. In fact, it can be used 
to determine the reflection factor f for the wall materials. 


FREE MOLECULE FLows AT SMALL Macuh NUMBERS 


Free molecule flow at low velocities has been studied 
for two cases: the flow around a sphere, and the flow in 
a tube. The calculation made by Epstein’ has shown 
that the drag coefficient Cp can be expressed 


Cp = 8/8V8/ry C(1/M) 


where C is a constant. Epstein determined C with 
diffuse reflection of molecules to be 1.442 for spheres 
that do not conduct heat and 1.393 for spheres that 
conduct heat perfectly. The experimental value of C 
for charged oil drop through air found by Millikan*! 
is 1.365. The lower experimental value can be ac- 
counted for by a small tendency toward specular re- 
flection of the molecules at the surface. By using the 
value of C = 1.393 and y = 1.405, the drag coefficients 
are indicated in Fig. 4 for small Reynolds Number Re. 
In this flow régime, the C, is independent of Re. Fig. 4 
shows that in the transition region between slip flow 
and free molecule flow, the drag coefficient actually 
decreases with decrease in Reynolds Number. 

The free molecule flow through a circular tube was 
first investigated by Knudsen'® and Smoluchowski.” 
Here the macroscopic velocity of the gas is constant 
across the section of the tube as could be expected from 
the large slip in the present case. The pressure drop 
can be expressed in terms of the pressure drop coefficient 
A defined by Eq. (30) as , 


A = (8V2/V2y)(f/2 — f)(1/M) 


(32) 


(33) 
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where f is fraction of diffuse reflection as defined pre- 
viously. By assuming f = 1, and y = 1.405, values of 
A are calculated from Eq. (33). They are indicated in 
Fig. 6 at the left side of the diagram. It is seen that, 
similar to the drag coefficient for the sphere, A actually 
decreases for a decrease of Reynolds Number in the 
transition region from slip flow to free molecule flow. 

The absence of Reynolds Number in Eqs. (32) and 
(33) simply means that the mean free path does not 
enter into this type of flow because of the smallness of 
the body dimension when compared with it. Thus, for 
both gas dynamical flow and for free molecuie flow, 
viscosity or Reynolds Number is a secondary effect. 
The primary effect is given by the Mach Number. In 
the slip flow region, both the Reynolds Number and the 
Mach Number are of the same importance. 


FREE MOLECULE FLow aT LARGE Macu NUMBERS 


Epstein and Smoluchowski assumed the velocity of 
macroscopic motion to be small in comparison with the 
microscopic molecular velocity in their investigations 
of free molecule flows. This means, of course, that the 
Mach Number of the flow is much smaller than unity. 
For large Mach Numbers, this restriction must be re- 
moved. The calculation should be simple because of 
the inherent simplicity of the physical situation. 
Sanger?’ 4 has solved this problem but his calcula- 
tion** is rather complicated. With the belief that a 
simpler derivation of the fundamental formulas will be 
helpful for an understanding of his results, the author 
will try to reproduce his basic equations by following 
Epstein’s general method for free molecule flows.’ 

If £’, n’, ¢’ are the components of velocity of a mole- 
cule in the directions x’, y’, z’ of a coordinate system in 
which the macroscopic velocity of the gas is zero, then 
since the velocity distribution is not modified by the 
impact with the body because of absence of collision 
between molecules in the free stream and re-emitted 
molecules from the surface of the body, the distribution 
is Maxwellian, or 

Nery? _ N(h/x)*/ 2 ge WG +0? + 5") (34) 
Nyzq’¢' is the number of molecules per unit volume in 
the range of velocities &’ to £’ + dé’, n’ to n’ + dn’, 
and ¢’ to ¢’ + df’ divided by dt’dn’dt’. WN is total of 
molecules per unit volume. h is related to the most 
probable velocity c; of the molecules in the free stream 
by the expression 


hk = l/c? (35) 
where ¢,2 = 2RT = 2(p/p). | 





Now, if the observer moves with a velocity ¢U, 
é.U, e3U in the directions x’, y’, 2’, then in the relative 
coordinate system x, y, 2 where the observer is con- 
sidered at rest, the velocities £, n, ¢ of the molecule in 


the directions x, y, z are 


—E= & —e¢U 
n=n7' — @U (36) 
f= - — eU 


Therefore, in the new coordinate system, the number 
of molecules having velocity component between £, 


n, § and + dt,» + dn, § + dé is 


h 3/2 
Ne,dé dn dg = N (“) em hE +eaU)*+ (nte2U)*+ & + e3U)2] 
U7] 
(37) 


Suppose there is a surface dS whose normal is the 
x-axis. During 1 sec., the molecules, having velocity 
components between &, n, f and & + dé, » + dn, § + df, 
and striking the surface dS, will be contained at a given 
moment in the cylinder with dS as base and the length 
V & + n? + £7 in the direction é, n, ¢. The volume of 
the cylinder is equal to the area of the base dS multiplied 
by the heights—£. The number of this kind of molecules 
is then N,,,dé dy d{(— &)dS. The number of molecules 
between &, n, ¢ and — + dé, n + dn, § + dé that will 
strike a unit area with x-axis as normal is then —£N,,,- 
dé dn df. 

To find the total number of molecules striking this 
unit area, one has to integrate —EN,,,d& dn df from 
&F=—o,7=—0,f=— wotot=0,7= ©, f= ow, 
The upper limit for ¢ is easily understood since no mole- 
cules with positive velocity in the x-direction can strike 


this area. Thus 


h 3/2 
ia=>- N (! 
Tv 
0 foo) @ . 
J af dn te —h[ (¢+e1U)2+ (n+e2U)?+(¢ +eaU)?] dt 


(38) 


The integration can easily be carried out; the result is 





1 U J 
n= WN a g hee + — {1 + erf(e,U vin| 
T 


(39) 
where erf(¢) is the error function™ defined as 
: 2 
erf(t) =—= e* ds (40) 
ax 0 


If the velocity U- of the motion of the area concerned is inclined with the angle @ with y-axis and in the x-y 


plane (Fig. 7), then 


é; = sin 6, 


é = cos 8, 


&3 = 0 (41) 


Hence 
with 


By m 
unit < 


Th 
e/E + 
be tra 
M.,'.¢ 


The r 
M.,’, 


bold 


For 
cules 
veloci 


By st 
presst 


Pi 
'/2 pl 





Eq. (4 


Fic. 7. 








SUPERAERODYNAMICS 661 


Hence, the number of molecules per second striking a unit area of the plate inclined at an angle @ to the stream 
with velocity U is 
1 ; U sin 6 ~ 
—h(U*) sin? @ PM ee 4 . 
=e + 1 + erf (U Vii sin 0 
av ah _ i ! ( ) 





n= 


By multiplying the above equation by the mass of the molecule, the mass m; of the stream per second striking a 
unit area of the plate is obtained, 


- { J2/¢42 sin? U . U.. 
mim Pee a 4 Vie sino 1+e(2 sno) (42) 


The component of velocity of the molecule in a direction, having directional cosines @’, é’, és’ with the axes is 
e;'—£ + e2'n + e3’¢. If this molecule is absorbed by the surface after striking it, the corresponding momentum will 
be transferred to the surface. For the surface with x-axis as normal as considered previously, the total momentum 


M.,’, 2’, ' pet second per unit area is 


3/2 0 es © ' 
Mea'e' = — (2) f dt an ff E(er’E + er'n + eg’te MEtAU + taly+ Gay) ge — (43) 


TT —@ 
The result of integration is 


M.,, e2’, es” 
é,’ 


1 , , , 
- 3 U? a (eres! + e262’ + eg¢s’)e~ 2” + E + €:(e1€1’ + e2€2’ + eres | [1 + erf (UVhe)}t (44) 








For calculating the pressure p, due to impact of mole- limiting casesof U = 0andU>c;, When U = 0, 


cules on a plate inclined at an angle @ to the stream of the impinging pressure is 

velocity U, 
Pi = "/4 pc? 

4 =sin@ ¢’ = —1 

é = cosd e' =0 (45) 


é; = 0 e;’ = 0 


This is one-half of the pressure acting on the surface 
calculated by the kinetic theory as the other half is the 
result of the reflection of molecules. For U > c, and 


By substituting Eq. (45) into Eq. (43), the impact 9 = 90°, the impinging pressure is 


pressure p; is calculated by 
Pi = p(?/ec? + U*) 





Pr = sj & — (U*/cé*) sin? 
> alaalh catia hil 
2 . \ 
(3 = + sin? ay] + erf (2 sin 0) (46) 


Eq. (46) reduces to particularly simple forms for the two 











p 4 
Fic. 7. Coordinate system for calculation free molecule flow. 


This checks with the result of Zahm! except for a factor 
of 2 which is accounted for by the fact that here only 
pressure due to impinging molecules is calculated, not 
the total pressure. Zahm obtained the result by simply 
assuming that U > c, and the molecules all having the 
same speed c, but with random distribution in their 
directions of motion. The later assumption is, of 
course, an oversimplification. However, if U > ¢, 
the spread of the speeds of the molecule is relatively 
unimportant. Then the present more accurate 
calculation should rightly give the same result as 
Zahm. 

To calculate the shearing stress r; due to impact of 
molecules on the plate, the direction cosines are 


@4=sind e' =0 
; 


€@=cosd e&' = —1 (47) 
é3 = 0 e;’ = 0 


By substituting Eq. (47) into Eq. (43), 
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cin wine = en (Ue) sin 0 
1/2 pU? U 


sin 6 cos 0 E + erf (2 sin 0) (48) 
L 

If the molecules reflect specularly from the surface, 
then motion of the impinging molecules normal to the 
surface is simply reversed in re-emission, while the 
tangential velocity is maintained. Then it is easy to 
see that the pressure due to re-emission is p, = p,; and 
the shearing stress due to re-emission is 7, = —7;. On 
the other hand, if the molecules re-emit diffusely, 7, = 
0, since no preferred direction exists. To calculate the 
pressure p, due to diffuse re-emission, one observes that 

according to the kinetic theory, 


p, = (mass emitted per second per unit area) (Va 2)c, 
(49) 
where c, is the most probable velocity of the molecule 
in thermal equilibrium at a temperature T, of the te- 
emitted gas. Thus 
c,? = 2QRT, (50) 
But the mass re-emitted must be equal to the mass m, 
striking the surface, as given by Eq. (40). Thus, 


Pre Vi Ge ON 
/ypU2 2 UU We 7 
. Vr G, U 
~~ 7 sin 6 E + erf (2 sin 6 | (51) 


Ci 





(diffuse reflection). 

The total pressure due to both impinging and re- 
emission, with f fraction of the molecules reflected dif- 
fusely and (1 — f) fraction of the molecule reflected 
specularly is 





= (52 


1/ .TT2 ” 2 _ ) x. ry? . 9 
iol? ~ @ ~ Dip, gue" 7,90 
The corresponding shearing stress is 

t/*/2 pU? = f(r:/1/2 pU) (53) 


These equations can be also written in terms of Mach 
Numbers by the substitution 


U ly 
a YM, = ly Yu (& ) 54 
Ci \; Cr " ) 


where J is the free stream Mach Number. 

The temperature 7, has to be determined by a con- 
sideration of the accommodation coefficient a. The 
energy content of the molecules from the wall is pro- 
portional (1 — f)7,° + fT, where 7;° is the stagnation 
temperature of the free stream. Therefore, according 
to the definition given in Eq. (20), 


Te? (A= fT° + fT) _ fTr°-T) 55) 
r=F Te-T, © 


where 7), is the temperature of the surface. If a and 


T,, are given, T, can be calculated from this equation. 
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Fic. 8. Surface temperature 7, for various angles @ to the 
wind in an atmosphere of 86 per cent Nz and 14 per cent O, at 
the density 1.9 X 10~*slugs percu.ft. Air temperature = 320°K. 
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Fic. 9. Impact pressure ratio p;/!/29V? for various angle @ 
to the wind in an atmosphere of 86 per cent Nz and 14 per cent 
O; at the density 1.9 X 10~* slugs per cu.ft. 


FREE MOLECULE FLOW OVER AN INCLINED PLATE 


To determine 7, for the case of an inclined plate 
moving through air, Sanger*‘ utilizes the condition of 
balance between energy lost through radiation and 
energy brought to the surface by impact of molecules. 
If E, is the energy radiated from the surface per second 
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per unit area, then it is related to the energy E; brought 
to the plate and the energy E, carried away from the 
plate by 

E,= E, — £, (56) 
Let €:ans» €or» €in, De’ the translational energy, the 
rotational energy, and the vibrational energy per until 
mass of the gas. Since experiments and theoretical 
considerations indicate that the re-emission is mostly 
diffuse, Sanger takes f = 1. Then, if m; is the mass 
brought to and re-emitted from the surface per unit per 
unit time, 


E; = m,(}, 2U? + €i trans. + trot. + Civin. 
E, Tr i(Erieans, ~ rot. a €rvib. 


(57) 
(58) 


The corresponding energies for a gas temperature 7), 
IS Evctrans.» €vror, ANA €y,;,- Sanger then assumes that the 
accommodation coefficients for translational and rota- 
tional energies are unity, while the accommodation 
coefficient for vibrational energy is zero because of its 
slow rate of adjustment. 


Then 
€itrans. ~ €r trans. trans. — €wtrans. 
trot. bie i €rrot. - €irot, na [wr ot. (59) 
Ewin. — Evin, = 9 
Hence, €rtrans, 4 €wtrans. €rrot. = €xrot. Ervin. = tvin. and 
c, = V2RT,, (60) 
Therefore, Eq. (58) can be written as 
E, = Mi Excans. + Cwrot. + Etvin.) (61) 


By substituting Eqs. (57) and (61) into Eq. (56), one 
has immediately 
E, Tt m,[*/2U? + ea. ee €wreans.) + 6... Exot.) ] 


For nitrogen and oxygen the equipartition value of the 
rotational energy is reached at even exceedingly low 
temperatures. Therefore, 


= 3/2RT 
2/2RT 


(62) 


€ => 


Hence Eq. (61) gives 
E, = m, [/2 U? + °/2R(T; — Tw)] (63) 


where m, is given by Eq. (42). 7; is the free stream 
temperature 7>. 

As an example, Sanger assumes that the air tempera- 
ture is equal to T; = JT) = 320°K. (47°C.). The air den- 
sity is 1.904 X 10~ slugs per cu.ft. The density ratio 
referred standard conditions is thus 0.802 xX 10-, 
corresponding to an altitude of approximately 55 miles. 
The mean free path according to Table 1 is thus approxi- 
mately 3.5in. Sanger actually assumed the air to be 
composed of 14 per cent oxygen and 86 per cent nitrogen 
without considering ionization and dissociation. Neg- 
lecting the radiation from the air to the surface, he as- 
sumed the emissivity of the surface to be 0.80. Thus 
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Fic. 10. Re-emission pressure ratio ~,/!/.»V? for various 
angle 6 to the wind in an atmosphere of 86 per cent Ne and 14 
per cent O, at the density 1.9 X 10~® slugs per cu.ft. 


E, = 0.800T p! 


where ¢ is the Stefan-Boltzmann constant for radiation. 
The results of calculation of the wall temperature 7, 
are shown in Fig. 8. It is interesting to note that be- 
cause of radiation, the wall temperature is considerably 
lower than the air temperature for low air velocities 
and suction side. The value of c, can be calculated by 
Eq. (60). The stresses acting on the surface can then be 
computed with Eqs. (52) and (53). The pressure ratios 
P:/1/2pU? and p,/!/29U? for impinging and diffuse re- 
emission (f = 1) are given in Figs. 9 and 10, also taken 
from Sanger’s paper. The shearing stress is given in 
Fig. 11. The extremely high shearing force is evident. 
If the dimension of a flat plate airfoil is much less 
than the mean free path, Eqs. (52) and (53) can be 
used to calculate the lift and drag. Thus, for diffuse 
reflection (f = 1), and for an angle of attack a, the lift 
coefficient C, and the drag coefficient Cp are 


C. = \2" 2{<) sin a cos a + 
y M\c; 


9 ct 
= 2 cos a erf (y2 M sin a) 


y M? ] (64) 


— 


_ — 

a | in? 2r 1 C ° 

c.. = }2 _* e— (7/2) M? sin a + pee (“) sin? a + 
. VY, M i M Ci 


| , oa 
(4, + 1) sin a erf (y? M sin «) (65) 





~ 


Sanger’s result for the flat plate airfoil is given in Fig. 
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+ Fro. 11. Shear stress ratio r/!/.9 V? for various angle @ to the 
wind in an atmosphere of 86 per cent Nz and 14 per cent O, at 
the density 1.9 X 10~* slugs per cu.ft. 


- - + « FF Pees 


12. The lift-drag ratio is less than unity, showing the 
poor efficiency of the surface under free molecule flow 
conditions. Of course, the calculated values hold true 
only for wing dimensions much less than the mean free 
path, i.e., less than 3.5 in. for the.case under considera- 
tion. 
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An- Investigation of the Airflow Underneath 
Helicopter Rotors 


ROBERT S. ROSS* 


Daniel Guggenheim Airship Institute 


SUMMARY 


The airflow underneath helicopter rotors, as determined by the 
classical blade element and momentum theories without the 
assumptions of constant chord or constant circulation along the 
span, was estimated from theoretically determined charts. 
These values were compared with those obtained from the actual 
measurements taken near a model rotor with a ‘‘hot wire” 
anemometer. 

Both time average and instantaneous readings of velocity 
were taken at simulated hovering and horizontal flight conditions. 
These measurements showed the theoretical estimation to be 
good at low angles of attack and at positions not too near the 
blade tips. They also indicated the pulsations that occur in 
horizontal flight, and showed these to be of unanticipated magni- 
tude. Positions of maximum and minimum effect were tested to 
get the range of pulsation. 

This method of determining rotor characteristics was different 
from those now employed, and the equipment used was built 
specially for these tests. The instruments were easy to handle 
and offered an improved picture of the distribution of rotor forces. 


DEFINITION OF TERMS 


= angle whose tangent is D/L 





a = 

a = angle of attack of the blade sections 

ay = angle between the rotor disc and horizontal plane 

B = angle the relative wind makes with the rotor plane 

b = number of rotor blades 

c = chord of rotor blade 

Cp = coefficient of drag of airfoil 

CL = coefficient of lift of airfoil 

D = drag of the rotor blade 

dD = drag for section of span dr 

aL = lift for section of span dr 

dr = elemental section of blade 

dT = elemental thrust for section of span dr 

dQ = elemental rotor torque for section of span dr 

dW, = elemental weight to be lifted by section of span dr 

K = drag over lift ratio 

Ki =1-—V1-— P? — 2KP 

L = lift of the rotor blade 

zx = {V1 — [(P? + 2sP + 2KP)/(1 — Ks)*)} 

M = (2s + P) 

n = r.p.s. of the rotor 

N = (1 — Ks) 

por P = ratio of the induced axial velocity to the rotational 
* velocity for hovering flight 

Q = rotor torque 

r = radius of blade to section 

R = tip radius of rotor blades 

p = mass density of air (0.002378 for standard air) 

s = ratio of the axial velocity to the rotational velocity 


Presented at the Rotating Wing Aircraft Session, Fourteenth 
Annual Meeting, I.A.S., New York, January 29-31, 1946. 
* Technical Director. 


C] = blade position with respect to the direction of the 
horizontal velocity 

u = rotational (or tangential) velocity of the rotor 

u’ = induced rotational (or tangential) velocity 

Ug = rotational velocity when rotor is tipped 

v = axial velocity of the rotor 

v’ = induced axial velocity 

v1’ = axial induced velocity for horizontal flight 

02" = induced axial velocity when rotor is tipped 

W = relative wind velocity for hovering flight 

Wo = weight to be lifted by rotor 

Y = horizontal velocity of helicopter 

Ye = horizontal velocity component in plane of rotor 
disc 

Z = resultant force on airfoil 

22 = resultant induced velocity from airfoil 

INTRODUCTION 


| Ee A GREAT NUMBER OF YEARS people have tried to 
solve the helicopter problems and build a practical 
man-carrying machine. However, it was not until the 
last couple of years that research has concentrated on 
the various helicopter difficulties and tried to find out 
more about them. Since the 15th century, adventurous 
people have been building helicopters of every variety 
and combination, each trying to make a completed 
machine. The basic principle of flight appeared so 
simple that no attempt was made to isolate any of the 
aerodynamics involved. As a matter of fact, few ex- 
perimenters realized that there were such problems as 
stability, vibration, or controllability. Their main 
trouble lay in trying to produce a rotor that would lift 
their machine with the power available.' 

Today, however, both rotor and propeller design have 
advanced to such a stage that the lifting problem has 
been solved. No sooner had a helicopter been built 
which would rise off the ground than all the complica- 
tions began to appear. It became evident that much 
concentrated research was needed on each particu- 
lar phase of helicopter flight to enable the design 
engineers to anticipate and try to solve these prob- 
lems. 

One of the unknown values which had come up time 
and time again was the distribution of the flow under 
the rotor blades. It was the purpose of this investiga- 
tion to try to give some answer to this question. In 
order to do this, it was necessary to develop experi- 
mental equipment capable of checking the validity of 
the use of basic helicopter theories. 
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DESCRIPTION OF APPARATUS 


All of the rotor tests were made in the 6 ft. 6 in. 
vertical wind tunnel of the Daniel Guggenheim Airship 
Institute. The rotor blades were mounted rigidly to 
the rotor shaft, the plane of the blades being 5 ft. 
above the tunnel mouth for hovering tests. The rotor 
shaft, in turn, was connected through a set of reduction 
gears to a '/3-hp., 10,000-r.p.m. 110-volt a.c. motor. 
The motor, gearing, shaft, and blades were all connected 
at the motor end to a specially devised arrangement of 
flat springs. The rotor end was supported by ball 
bearings mounted inside of a protective tube that ex- 
tended from the blades to the motor housing. 


Force Apparatus 


The entire arrangement, except for the rotor blades, 
of course, was housed in a streamlined shield to prevent 
any interference from the air and to reduce to a mini- 
mum any effect of the machine on the airflow itself. 

The motor was prevented from moving in either an 
axial or rotating direction by the flat spring support on 
the bottom of the instrument. Two of the springs were 
crossed; this enabled the motor to rotate a fraction of a 
degree but prevented it from moving in the axial direc- 
tion. The other spring prevented the first set of springs 
from rotating but allowed them to move axially just a 
fraction of aninch. By mounting strain gages on these 
two sets of springs, the tendency to move in either the 
axial or rotary direction could be measured. These 
values, when calibrated, gave the thrust and torque of 
the rotor blades. 

The strain gages, which were mounted on the springs 
in this force apparatus, formed part of a bridge circuit. 
The unbalance of this circuit due to the stretching of the 
gages was amplified and recorded on Esterline Angus re- 
corders, which were calibrated to measure the thrust 
and torque of the rotor blades. This part of the ap- 
paratus therefore gave the performance characteristics 
as measured at the shaft and was used as a qualitative 
check on the induced velocity measurement. 


‘“*Hot Wire’’ Indicators 


In order to measure the velocity and direction of the 
air under the rotor blades, a special ‘‘hot wire’’ instru- 
ment was used. These wires were mounted at the end 
of a long probe, which, in turn, was supported in an 
adjustable clamp enabling its placement at any desired 
position under the blades. The ‘“‘wires’’ themselves 
consisted of two separate units: one for the direction 
indication, and the other for the velocity. The direc- 
tion unit was made up of two '/s-in. long wires in the 
form of a V. The velocity unit was one !/s-in. long 
wire set at right angles to the plane of the direction V. 
All wires used were 0.0006-in. diameter platinum 
strands. They were strung between needle points 
properly spaced at the end of the adjustable probe. 


The two sides of the V of the direction wire made up 
two arms of a bridge circuit. If the bridge showed a 
balance, it meant that these two sides of the V were 
equally cooled by the airflow, and the direction of this 
flow could be read on a large protractor on the probe. 
The velocity wire made up part of another bridge cir- 
cuit, the unbalance of which was amplified and recorded 
on another Esterline Angus recorder. This recorder 
was calibrated to read the velocity of the air passing 
over the wire. The entire “hot wire’ unit was no 
larger than a sphere of approximately */,-in. diameter. 
This enabled ‘‘point’” investigations, since it was so 
small. 

The “‘hot wire”’ instrument was calibrated in a special 
small wind tunnel. It consisted of a 7-in. diameter 
tube about 3 ft. long, at one end of which was mounted 
the motor and fan. The tube was fitted with screens 
at the entrance and “‘egg crates’’ at the exit so as to ob- 
tain relatively smooth flow. This small wind tunnel 
was equipped with a 1/j-in. diameter pitot tube con- 
nected to an alcohol manometer inclined to 5° with the 
horizontal. A small hole in the side of the tunnel !/, in. 
in diameter enabled one to insert the ‘‘hot wire’ probe 
until the wires rested at the center of the tube. Veloci- 
ties from 0 to about 30 ft. per sec. could be obtained 
here. The speed of the tunnel motor was regulated with 
a variable inductance. Both the velocity and direction 
wires were calibrated in this tunnel for low velocity 
readings. 

In order to measure the instantaneous velocities, a 
“hot wire” about '/2 in. in length was installed at the 
end of the probe. This wire was connected through a 
bridge circuit to an oscillograph on which a photo- 
graphic record of the velocity changes could be made. 
A two-element unit was used so that a 60-cycle timing 
wave could be photographed alongside the velocity pic- 
ture. Records were made by either drawing the film 
over the light spots in an X direction while spots were 
restricted to move in the Y direction or by allowing the 
light spots to scan a photographic film screen once 
while moving in both the X and Y directions. The 
film, when developed, made a permanent record that 
could be easily studied. 


Types of Rotor Blades 


The rotor blades used for the tests were made of 
mahogany wood shaped to an NACA 0015 profile. 
Three sets of rotors were tested, each set employing 
three 2ft blades. The first set had a constant chord of 
2 in. and was twisted so that the root was ‘10° greater 
than the tip. The second and third sets were tapered 
so that the root chord length was 2°/, in. and the tip 
was 1'/, in. One of these sets was not twisted, while 
the other had a:10° twist similar to that of the constant 
chord set. The blades, which were fitted with a tapped 
socket at the root end, were screwed to a small hub and 
held rigidly in place by a lock nut. 
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TEST PROCEDURE 


‘*Hot Wire’’ Measurements 


The first investigations made were those with the 
“hot wire.’’ Before starting the calibrations, it was 
necessary to check the ‘‘hot wires’’ themselves to make 
sure that they were well soldered to their needle sup- 
ports and that good electrical contact was made 
throughout the circuits. The “hot wire’ unit was then 
plugged into the long probe, and the probe was swung 
around into position so that the wire rested in the cen- 
ter of the small calibrating tunnel. The electrical ap- 
paratus was then calibrated so that the Esterline Angus 
recorder would read the velocity over the wire, while 
the angle of the direction wire was maintained at zero 
by rotating the entire probe. In order to best utilize 
the recorder for the majority of the velocities to be 
encountered, it was necessary to balance the velocity 
bridge with the tunnel operating at a low velocity. 
This gave better accuracy to the values taken but pro- 
hibited the recording of zero velocity itself. The ‘‘hot 
wire’ was calibrated before each run to make sure that 
such things as the change in temperature of the air in 
the room and its effect on the wire did not enter into the 
results. 

So as to produce hovering flight conditions, the rotor 
blades were placed in the tunnel so that the plane of the 
blades was horizontal and about 5 ft. above the tunnel 
mouth. Because it was a vertical tunnel, the opening 
under the blades helped remove the ground effect. The 
turning vanes and corner of the tunnel, the first solid 
obstructions under the blades, were about 25 ft. from 
the helicopter rotor. All of the tests were run in semi- 
darkness. A Strobotac was used to determine the rotor 
speed, which was kept constant at 450 r.p.m. This 
value was chosen arbitrarily because the equipment ran 
well at that speed. 

After the “hot wire”’ was calibrated, the tunnel room 
was darkened, and the probe was swung into a position 
so that the ‘“‘hot wire’’ could be easily moved under- 
neath the rotor disc. Velocity and direction readings 
were then taken in a plane about 2 in. below the rotor 
disc, and a complete traverse was made from a radius 
of about 2 in. to about 24in. At extremely low veloci- 
ties the readings could not be taken, and at positions 
where the air was turbulent the readings were not ac- 
curate. However, regions of these conditions did not 
contribute much to the effectiveness of the blade, so 
this was no great loss. Readings were taken radially 
at l-in. intervals, the resultant velocities and direc- 
tions being recorded. 


Horizontal Flight Tests 


So as to obtain horizontal flight conditions, the rotor 
was inclined in the air stream until the plane of the 
blade was parallel to the direction of the wind-tunnel 
air flow. The supports of the motor housing were 


moved back so that the rotor blades themselves re- 
mained well within the air stream. The rotor chosen 
to be investigated was the one with the nontapered, 
twisted blades that were set at (8 + ay) = 20° at the 
root. These were mounted rigidly at the rotor hub, and 
the “hot wire’ probe was mounted so that measure- 
ments could be taken at 2 in. below the plane of the 
blades. The velocity and direction ‘‘hot wire’ unit 
was first used, and readings were taken at 1-in. intervals © 
along the radius of the blades. Then the instantaneous 
velocity wire was installed and readings were taken only 
at several representative points. 

The wind-tunnel velocity, which was determined by a 
pitot tube mounted in the air stream, was held at a con- 
stant value depending. upon the test desired. Since 
the wire would burn out at zero air velocity, care had 
to be exercised that air was flowing past the wire before 
the current was turned on. This was necessary because 
high current had to be maintained through the wire to 
obtain the sensitivity desired. 


HELICOPTER THEORY 


An examination of the airflow around a helicopter 
rotor blade with the primary purpose of determining 
design criteria was the main function of this investiga- 
tion. By the use of the well-known blade element and 
momentum theories, formulas were developed which 
would give the characteristics of the rotor blades under 
various attitudes of flight. These formulas, which be- 
come quite complicated when both the lift and drag of 
the blade element are considered, were plotted in curve 
form and in that way present a simple and rapid 
means of their solution for any particular flight condi- 
tion. 

All the important design factors were included in one 
“design coefficient’”’ (C,bc/r) which is dependent upon 
the spanwise blade loading and the rotational speed of 
the blade. The physical characteristics of the blade 
shape would therefore be determined by the structural 
ability of the blade to handle the load given it. Once 
the possible blade loading is decided upon, it would be 
a matter of applying these values to the above- 
mentioned curves to obtain the C,bc/r value for any 
section. 

The general solution of the problem was considered 
first, and then various possible flight conditions, such 
as hovering and horizontal flight, were investigated. 
The entire theory was determined primarily from the 
basic blade element and momentum theories and was 
developed from there. 

Taking a rotor blade section of span dr (as illustrated 
in Fig. 1) and examining the airflow at that point where 
r is the radius of the rotor to that section, ¢ is the chord 
of the rotor section, and R is the tip radius, one finds the 
conditions illustrated in Fig. 2. These can be analyzed 
as follows: 
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Fic. 1. Rotor blade definition. 
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Fic. 2. Velocity and force diagram on cross section of rotor 
blade. 





Let the terms be defined so that 


v = the velocity entering the rotor disc (axial 
velocity of the rotor) 

u = the rotational velocity of the rotor 

v’ = the increase in axial velocity or induced 
axial velocity 

u’ = the increase in tangential (rotational) ve- 


locity or induced tangential velocity 
dL = lift of the rotor section 
dD = drag of the rotor section 
W = velocity relative to the rotor section 


If W is the relative velocity of the air with respect to 
the blade section, then the lift component dZ must be 
perpendicular to W and the drag component dD must 
be parallel to it. The resultant of these lift and drag 
components is Z, which makes an angle a with the lift 
vector. 

This force Z has a reaction on the air equal and op- 
posite to it. Call the velocity set up by this force an 
induced velocity and give it the value Z,. Z: can then 
be broken up into its horizontal and vertical compo- 
nents u’ and v’, respectively. If 8 is the angle the rela- 
tive velocity W makes with the horizontal axis, then 
Z, makes the angle (8 + a) with the vertical axis, and 


tan (8 + a) = u’/v’ (1) 


According to the momentum theory, if the overall 
increase of velocity is v’, then the increase at the blade 
is one-half of that or v’/2. Applying this to both the 
tangential and axial velocities, it is found that the rela- 
tive velocity W is made up of the axial velocity 
[v + (v’/2)] and the tangential velocity [u — (u’/2)]. 
From this new vector triangle, it is clear that 


tan 8B = [v + (v’/2)]/[u — (w’/2)] (2) 
Expanding Eq. (1), 
u'/v’ = (tan 8+ tan a)/(1 — tan Btana) (3) 


Since the tan a = dD/dL or Cp/C,, Eq. (2) can be sub- 
stituted in Eq. (8) to give 


uv + (0’/2)]/[u — (u’/2))} % (22) (4) 


v1 — {fv + (0’/2)]/[u — (u’/2)]} | \G/Gr 


Let Cp/C, = K and simplify Eq. (4). This gives 


,| (20 + v’ 15 gig jo Saree 
’ [G@ts)+«]-« E sages 9 ©) 


Expanding Eq. (5) gives 








Quv’ + v’? + 2Kv’u — Ko'u’ = Quu’ — u’? — 
2vu’'K — u’'v’K (6) 
Collecting terms and rearranging Eq. (6) so that it 
would be possible to solve for u’ yields 
u’? — 2u’(u — vK) + v’? + Qo’ + 2Kv’u = 0 
(7) 
Solving Eq. (7) for uv’, 
u’ = 2(u — vK) = 
V4(u — vK)? — 4(v’2 + ov’ + 2Kv’u) (8) 
u’ = (u — vK) + 
V(u — vK)? —(v'? + 2Qvv’ + 2Kv'u) (9) 
In order to analyze Eq. (9) to find out which sign ap- 
plies, rearrange it into this form: 
u’ = (u — vK) + (u — vK) X 
ne Pe’ ons 12 , / 
yi vK) (v’? + 2vv’ + 2Kv’u) (10) 




















(u — vK)? 
° ai : v’? + Qu’ + 2Kv'u | 
u’ = (u — vK) E “1 = Gu 0K)? | (11) 


To simplify the equation even further, consider the 
condition where the axial velocity v is equal to zero. 
The expression (u — vK) then becomes u, and the 
value inside the brackets must be equal to or less than 1 
because u’ cannot be larger than u. This means that 
the negative sign must be employed, since the value 
under the radical always lies between zero and 1. 

Consider now the force that can be obtained from this 
blade section. Since ; 
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force = thrust = mass X acceleration = 
(mass/sec.) (change in velocity) 


the following relationship can be set up: 
dT = 2arpdrv average axial) x U (axial) (12) 


where dT is the elemental thrust force and p is the mass 
density of the air. The average axial velocity is equal 
to the velocity at the section or [v + (v’/2)], while the 
change in axial velocity is the v’ mentioned earlier. 
Substituting these values in Eq. (12) gives 
dT = 2xrdrp[v + (v’/2)]v’ (13) 
which, when simplified, yields 
dT /dr = xrp(2vv’ + v’?) (14) 
To find the axial thrust of the rotor section by 


another method, the following equation can be set up: 


thrust = L cos B — Dsin B (15) 


or 
! velocity)? 
thrust = (C, cos B — Cp sin B)p ( a area (16) 


where the velocity equals the relative velocity (W) and 
the area equals the area of the rotor blades. Eq. (16) 
can then be written for the elemental sector considered 
previously: 
dT = (C, cos B — Cpsin B)(p/2)W? bedr (17) 
where } equals the number of rotor blades used. From 
Fig. 2 it can be seen that 
W = [v + (v’/2)]/ sin B (18) 


Substituting this value in Eq. (17) and rearranging the 
equation to solve for d7‘/dr results in 


dT = . ov (P\| v + i 
oe (Cz, cos B Cp sin B) (2) = bc 
(19) 


Eqs. (19) and (14) both give values for d7/dr and can 
therefore be equated: 


2rr(2vv’ + v’*) = (C, cos B — 


Cp sin 8) | 


v + (v’/2) 
sin B 


2% 
[ve (20) 


Since Cp/C, = K or Cp = KC,, Eq. (20) can be rewrit- 
ten so that 


2rr(2vv’ + v’?) = C,(cos B — 


K sin B) pte”) al be (21) 


sin 8 
Solving Eq. (21), 


Crbc = 2r(2vv’ + v’?) 
r (cot 8B — K)[(2v + v’)2/4 sin B] 


This can then be put into the simpler form: 


(22) 








CLbc Sx sin 8 


= - ———__——- (23) 
r (cot 6 — K)[2(v/v’) + 1] 








which is the final design formula. 


Condition of Hovering Flight 

Now that the design Eq. (23) is available, apply it 
to the simplest condition, hovering flight, where the 
axial velocity of the rotor blade is v = 0: The rota- 
tional induced velocity u’ from Eq. (11) becomes 





u’ = u{l — V1 = [(v’? + 2Kv'u) u?}} (24) 
and the tan 8 from Eq. (2) becomes 


v’/2 


tan B = on 
Uu vy’? a 2Kv'u (25) 
S—-—j | = Os es. 
2 u? 


The new equation for the relative velocity W becomes 


W = y(5) 4 E as (1 i y' fy id 
2 2 u2 


(26) 














from the original equation taken from Fig. 2, where 





W = Vivo + (v’/2)]? + [u — (u’/2)]? (27) 
The sin 8, which hecomes 


sin B = 





fp? [uf | v®4+9Keu\} (28) 
ti -4- ==) 


can be simplified to the following form: 





m / eT au one ge ey EY = - 
sin B = v’/\2u? + 2u(V u? — v’? — 2Ko'u — Kv’) 
(29) 





Now that v = 0, the new design formula, Eq. (23), be- 
comes 


C,bc/r = 8m sin B/(cot B — K) (30) 


Substituting both Eqs. (29) and (25) in Eq. (30) results 
in 














Crbc _ se| Tina + Vu? — 9’? — 2Kv'u - =a 
r [(u + Vu? — ov? — 2Kv'u)/v'] — K 








(31) 
This equation can be simplified by collecting terms to 


C,bc Sav’? 








(32) 


By dividing the numerator and denominator by u’*, 
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Sr(v’ /u)? 


Vf2 [1 — K(v’/u) + V1 — 2K(v’/u) — v 


9 


2, u?| 


Since the ratio, v’/u appears so often, let it be called P, 


and the design Eq. (33) for hovering flight becomes 
Cybe _ IPE i) to 
r (1— KP + V1 — 2KP — P?*)” 


All the equations have now been derived whereby one 
can obtain the design condition for hovering flight. If 


the weight to be lifted is Wo, then a dW can be assigned 
to each dr of the blade and replaced for dT in Eq. (14). 
If v is zero, as it is in hovering flight, then this equation 
becomes 

dW,/dr = mrpv’? (35) 
and if a value is assigned to dW,/dr, v’ can be determined 
at any given radius. If the distribution of the in- 
duced velocity is uniform over the rotor disc, Eq. (35) 
can be written in the integral form and solved: 


Wo = mov’? Jo* rdr (36) 
Wo = mpv’?(R2/2), (37) 
or 
vo’ = V2Wo/mpR? (38) 
v’ = (1/R)V2Wo/xp (39) 


This equation gives an average v’ for the entire disc. 

However, it is impossible to make v’ constant through- 
out the disc, since, if it were, the value P, (v’/u), 
would become extremely large near the center of the 
rotor. By examination of Eq. (34), it is clear that P 
can never be over 1, since, if it were, the value under 
the radical sign would be negative and no design coeffi- 
cient C,bc/r could be obtained. It is therefore apparent 
that in the actual design v’ will probably decrease as r 
approaches the center of the rotor. 

If a rotor loading is assumed in which the load is uni- 
form throughout the disc—that is, that each dr carries 
the same weight—then the induced velocity for any 
section dr between R; and R; would be 


Wo = In mrpv’? dr (40) 


ite ele =] (41) 
. Tp V Ry? _— R? 


Since V 2Wo/ mp is a constant, v’ just varies with the 
reciprocal of VR? — R,*. As the values of R. and R, 
approach the center, this radical would decrease in 
value, while the values for v’ would again increase. 
This also proves to be the wrong system of loading, 
since no blade could be designed to furnish the rela- 


tively high values of v’ occurring. 


or 


(34) 


Assume then a v’ that is proportional to 7, say fr. 
If this value is now substituted for v’ in the equation 
pb = v’/u and if uw is given the value 2zrn, then 

P = fr/2nrn or f/24n (42) 


which is a constant for any r. Substituting this new 
value for P in Eg. (34) gives C,bc/r as a constant also. 
This means that if the C, does not vary through the 
spare of the rotor blade, then the chord C of the blade 
will increase with the radius r. This would indicate a 
blade narrow at the root and wide at the tip. Aero- 
dynamically, the blade design is now possible, but, 
structurally, it does not appear to be plausible. Either 
a different induced velocity distribution is necessary, 
or the C,; value will also have to change along the span. 


Use of Charts 


So as to determine the theoretical velocity distribu- 
tion underneath the rotor blade in hovering flight, three 
curves, Figs. 3, 4, and 5, are used. They are, namely, 
a graph of the rotor blade characteristics; a graph of 
the value C,bc/r vs. P; and a graph of 8 vs. P. To 
illustrate the use of these graphs, a particular case is 
chosen and the theoretical distribution is determined. 

Assume a 4-ft. diameter rotor with three blades of 
2-in. constant chord, each rotating at 450 r.p.m. Let 
the blade itself have a 10° twist, the tip section being 
set at 10° with respect to the horizontal plane and the 
root being 20°. Use as blade section an NACA 0015 
Airfoil. Sections along the span 4 in. apart will be 
chosen so that those to be considered are at r = 0, 4, 
8, 12, 16, 20, and 24 in., respectively. Choosing the 
section at r = 12 in. as an example, the procedure for 
determining the values v’ and u’ at that point are as 
follows: 

Since it is not known just what C, exists at r = 12in., 
three assumptions could be made and the actual point 
found by plotting these values on Fig. 6. Values of 
C, = 0.85, 0.51, and 0.27, were arbitrarily chosen. If 
C, were 0.51 from Fig. 3, the angle a) would have to be 
7° and the L/D would have to be 20. Since the blade 
has a 10° twist and is set at 20° at the root, the angle 
(8 + a) at r = 12 in. would be 15°. This would 
therefore mean that 8 must be 8°. From Fig. 5, if 
B were 8° and the L/D ratio were 20, then P would 
necessarily be 0.272. However, the C,bc/r value with 
C, as 0.51 at a radius ry = 12 is equal to 0.255. From 
Fig. 4 this gave P equal to 0.199. These points were 
then plotted on Fig. 6, and the same procedure was fol- 
lowed through for C, = 0.85 and 0.27. By this 
method, it was determined that P = 0.218 atr = 12. 
If the rotor is operating at 450 r.p.m., then the tangen- 
tial velocity (u) of the blade section considered is equal 
to 2mrn or 47.1 ft. per sec. P was defined as v’/u, s0 


v’ is equal to Pu or 0.218 X 47.1 = 10.27 ft. per sec. © 


This value was plotted on Fig. 8 as the value of the axial 
. . . / 

velocity. In a similar manner the other values of ? 

were determined and the entire curve was plotted. 
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Fic. 3. NACA 0015 airfoil characteristics. 
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Incidence chart for various drag/lift ratios. 
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Fic. 6. Graphic determination of velocity ratio. 
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Induced velocity distributions for twisted constant chord 
rotor at 15° (hovering flight). 

















Induced velocity distributions for twisted constant chord 
rotor at 20° (hovering flight). 
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Fic.9. Induced velocity distributions for twisted constant chord 
rotor at 25° (hovering flight). 
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Fic. 10. Induced velocity distributions for twisted constant 
chord rotor at 30° (hovering flight). 
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Fic. 11. Induced velocity distributions for nontwisted, tapered- 
chord rotors. 
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Fic. 12. Induced velocity distributions for twisted, tapered- 
chord rotors. 


To determine the magnitude of u’ the value of P = 
0.218 was used in the formula u’ = u[1 — Vi-~- Fa 


KP]. K was determined by finding the C, of the 
section equal to 0.62 from Fig. 6 and the K value at 
this condition equal to 0.0532 from Fig. 3. This gave 
u’ = 1.75 ft. per sec. 

In order to show that these curves in Figs. 3, 4, and 
5 could be used for the determination of v’ and w’ for 
blades at other conditions and of other shapes, this pro- 
cedure was continued for this same blade at (8 + ao) = 
30°, 25°, 15°, and for both tapered untwisted and 
tapered twisted blades at (8 + ao) = 25°, 20°, 15°, and 
10°. These are shown in Figs. 7 through 12, respec- 
tively. 





Condition of Horizontal Flight with v = 0 


At the extreme condition of flight, where the heli- 
copter is moving forward and the plane of the rotor re- 
mains horizontal, the value of v’ and uw’ varies with the 
blade position as well as along the radius.* However, 
the same values of P as those in hovering flight are 
still valid, except that they are now defined differently. 
If the forward velocity of the helicopter is Y ft. per sec. 
then P becomes equal to v’/(u — Y sin 6), instead of 
v’/u as it had been in hovering flight. 


If the helicopter is considered as traveling forward’ 


with a velocity Y and yet the plane of the blades re- 
mains horizontal, then the new rotational velocity at 
any particular radius would be (uw — Y sin @), where é 
is defined as illustrated in Fig. 13. 
Since there is no vertical velocity, this condition 
represents a helicopter in hovering flight but moving 
horizontally. (There is no change in altitude.) The 
equation for this condition is still Eq. (34), except that 
p would now be defined as 1;’/(u — Y sin @), where 0 
is the new induced axial velocity. This axial velocity 
is now a function of the blade position @ and will vary 
from a maximum at 0 = 3x/2 to a minimum at 6 = 2/2. 
This means that if p is found by means of Figs. 3, 4, 
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Fic. 13. Top view of rotor disc in horizontal hovering flight. 


and 5 for the hovering condition, then this value of p 
for any particular radius r could be used in the deter- 
mination of the new 2’: 


v1’ = p(u — Ysin 6) (43) 


Considering now a small section of rotor span dr 
and solving for the lift produced by it, it is found that 


dL/dr = prd6(v;'?/2) (44) 


This comes from the momentum formula, which states 
that 


lift = (mass/sec.)(change in velocity) 


If it is desired to find the entire lift produced by this 
section dr at radius 7, then 





dL : ” ” 
a iF i (u — Y sin 0)*d6 (45) 
r 2 Jo 
or 
dL pp {7 , oe 
= on f (u2 — 2uY sin @ + Y? sin? @)d6@ (46) 


This gives 


dL op? @ _ sin 20) )* 
ai E + 2uY cos @ + Y? (5 - =) (47) 
0 


dr 2 2 


dL/dr = (p*rp/2)(2ru? + Y*r) (48) 


If the helicopter were not moving horizontally, then 
this equation would read 
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Fic. 14. Side view of inclined rotor disc in horizontal flight. 





dL/dr = (p*rp/2)2ru? (49) 


This formula therefore shows an increase of lift for the 
section considered equal to (p*rp/2)(Y*r). This in- 
crease of lift, of course, was based on the assumption 
that the rotor blades would continue to rotate at the 
same speed as they had been in hovering flight. The 
amount of torque required for a similar section should 
therefore be examined for both hovering and horizontal 
flight, as just described, to determine what increase of 
power would be required to maintain this rotational 
speed, 

The formula for the torque required‘ for section dr 
at radius 7 is similar to the thrust formula, Eq. (48), 
except that the value p is replaced by Ki, which equals 


(1 — v1 = bp? — 2Kp), and a distance factor (r) is 





added. 
dQ/dr = (Kippr?/2)(2xu? + Y*r) (50) 


represents the condition of hovering forward flight, 
while 


dQ/dr = (Kippr?/2)2xu? (51) 


represents the condition for still hovering flight. 
(Kippr?/2) Y*x, the difference between these two values, 
is the increase of torque required for a horizontal move- 
ment of the blades at a velocity Y over the torque re- 
quired for plain hovering flight for the same rotor 
speed. Therefore, although a gain of lift of (p?pr/2) X 
(Y*r) is felt during this horizontal movement, an in- 
crease of torque of Kippr?(Y*r)/2 was necessary in 
order to produce this result. 


In normal horizontal flight the axis of the helicopter rotor is usually tilted in the direction of motion. This 
causes a component of the total lifting force of the rotor to be directed horizontally and enables horizontal motion 
without the need of auxiliary rotors. This type of flight condition can be investigated with the formulas just 
completed. From Fig. 14 it can be seen that velocity v now equals Y sin a and is the same throughout the disc. 


The new rotational velocity is equal to uw. = Y2 sin 6, where Y2 = Y cos ay. 


Therefore, 





“ 
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Uy = u — Y cos q sin 8 (52) 


The thrust of the rotor then can be written for a section dr at distance r: 


dT/dr = pr Dg [Y sin a, + (v2’/2)]v2’ dé (53) 
where 
vo’ = P(u — Y cos qy sin @) (54) 
1T ee , f , : 
= pr J {Y sin a + [P(u — Y cos q sin 6)/2]}[P(uY<cos a sin 6)]dé (55) 
dr 0 
or 
dT/dr = prPr[2Yu sin a, + (P/2)(2u? + Y? cos? ay) 56) 
In order to check this equation, let Y = 0; then 
dT /dr = prP[{ + (P/2)(2ru?)] = (prP?/2)(2ru?) 57) 


which is the same as Eq. (49). Then let a; = 0 and 
dT/dr = (prP?/2)(2ru? + Yr) (58) 
which is the same as Eq. (48). 
The lift of the rotor blade for the section considered now decreases according to the angle a. 


dL/dr = (dT/dr) cos ay (59) 


The remainder, (d7‘/dr) sin a, goes into the forward thrust of the rotors. The forward thrust equation for the 
rotor at section (r, dr) is 

(dT /dr) sin a, = prPxr sin a[2VYu sin a, + (P/2) (2u? + Y? cos? a)] (60) 

In this equation all the values are known except P, and that can be obtained from the charts on either hovering 

or vertical flight conditions. These values may then be summed up for all the dr’s and the total forward thrust 


obtained. 


Condition of Vertical Flight 

Previously, the specific case of hovering flight was considered, but the formulas developed for that condition did 
not take into account what would happen in the event the ship would rise or descend. Eq. (23) is a general equa- 
tion for any condition and can be solved provided the terms in it can be determined. This was done in Eq. (34) 
for hovering flight but becomes a little more complicated for vertical flight. 

If the vertical velocity is equal to v as defined in Fig. 3, then, 


i. ws +! /9 o o//9* 
' v + (v’/2) (v + v’/2) 
sin B = : geo (61) 


The term wu’ can be defined by Eq. (11) and substituted in Eq. (61) to give 


l 








“Sry ‘ ra >) 9 (R9 
v’? + Qvv’ + 2Kv =\T (62) 


| 
1 + “ 2u.— (u — vK) (: — @/1 —- 7 — 
(2v + = (u — vK)? 


In the same manner, the cot 8 cam be found to be 


\ u—vkK gf? + 2vv’ + 2Kv'u 
cot B = Ng ( ; ) 1- l— - a rar , (vy + (v’/2)] 63) 
f (u — vK)? 


These terms, as defined in Eqs. (62) and (63) can be substituted in the general Eq. (23) to give the solution de- 
sired. Through simplification, this equation can be put into another form. By dividing both the numerator and 
denominator by uv? and letting v’/u = Pandv/u = s, the equation becomes 


SrPM 


C,be 
: = (64) 





L,)]? (NQ + Li). — KP] 


r VM?+(2- M1 — 





where M = (25+ P),N = (1—Ks) andl; = V1 — [(P? + 2sP + 2KP)/(1 — Ks)?]. 


This Eq. (64) is now in a form that is general enough to be applied easily to either hovering, vertical, or hori- 


zontal flight. 
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RESULTS 


Hovering Flight 


By employing Figs. 3, 4, and 5, the theoretical in- 
duced velocity distribution of three types of blades at 
various angles of incidence were calculated and plotted 
by the means described in the section on the ‘‘Use of 
Charts.’ The results of these calculations were plotted 
in Figs. 7 through 12 and show the relative values that 
can be expected when the various shapes and settings 
are used. 

As would be expected, the curves show that there is a 
gradual decrease in the induced velocity as the blade 
angle of incidence is decreased. This amount is most 
easily observed in Figs. 11 and 12, where the group of 
blades are plotted together. By comparing the curves 
of the 20° blade settings, it appears that the twisted 
blades have the more even distribution. The tapered 
blade, which is not twisted, carries an even greater load 
toward the tip than the constant chord blade. The 
two twisted blades have a similar distribution, the 
tapered blade being a little more evenly loaded than 
the other. This tends to indicate that the tapered 
twisted blade would probably be the best of the group 
structurally, while the constant chord blades, being 
simpler to construct, are a close approximation. 

The induced velocities of the constant chord rotors 
were measured in the wind tunnel by means of the “‘hot 
wire’ units discussed earlier, and the experimental 
points obtained were plotted on Figs. 7-10. The ex- 
perimental values for the 15° and 20° blade settings are 
fairly close to those predicted by the theoretical curves. 
However, as the angle of incidence of the blade is in- 
creased, the relationship between experimental and 
predicted values becomes less apparent, especially to- 
ward the tip. Since no tip correction was incorporated 
in the theoretical curves, it was to be expected that a 
discrepancy would appear. However, the breakdown 
of the induced velocity distribution does not localize 
itself to the tip but starts about at the center of the 
blade. By comparing the four constant chord blade 
curves, it appears that this breakdown acts almost as 
an ordinary stall but settles down to practically the 
same values for the three highest angle settings. 

“The measured rotational induced velocities were 
fairly close to the predicted values in’ those regions 
where the axial induced velocities checked. Near the 
tip, however, where the measured axial velocities be- 
gan to decrease, the rotational measurements showed 
a slight increase. This indicates that the air has be- 
gun to follow the blade around as if it were a paddle in- 
stead of a propeller and is an evidence of wasted power. 


Accuracy of Measurements 


The velocities close to the hub were not taken be- 
cause the equipment was set to measure only velocities 
above a certain base value. As was explained earlier 


in the discussion on the apparatus, this enabled more 
accurate readings of the velocities that were taken. 
The tip velocities also tended to zero and were omitted 
either because they were too low or because the tur- 
bulence in that region was too great to be assigned a 
constant value. The remainder of the velocities were 
measured to within an accuracy of about | per cent. 
This figure is an estimate, since it depends mainly on 
the ability to read the record and get a good average 
value for the velocity indicated. The actual velocity 
records, however, oscillated some 2 mm. out of 1S0, 
necessitating the use of an average value. The wind- 
tunnel air, when not disturbed by the rotor blades, 
would record as a straight line, the fluctuations in it 
being too fast and too small for the Esterline Angus 
meter to follow. The higher velocities under the rotor, 
therefore, could easily be measured with a better than 
1 per cent accuracy. 

The angle measurement, which was also taken, pre- 
supposes that there is not much radial velocity present 
and that the majority of the air is traveling in a plane 
perpendicular to the radius. The V-wire, therefore, 
measures the angle in that plane, and the velocity wire 
measures the velocity in that plane. On the whole, 
this assumption of low radial velocities is valid for the 
rotational velocity used and was checked by making a 
survey with the “hot wire’’ probe along that radius 
which would cause the V-wire to lie in a plane made up 
of the axis of the rotor and the radius. 


Horizontal Flight 


One rotor arrangement was chosen to be investigated 
thoroughly, and it consisted of the three constant chord 
blades set at 20° with respect to the rotor disc at the 
hub. With these blades revolving at 450 r.p.m. and 
inclined in the wind tunnel so that the plane of the disc 
was parallel to the air stream, several surveys were 
made under the rotor at various wind-tunnel velocities. 
Both average velocity surveys and instantaneous rec- 
ords were taken at positions that would be equivalent 
to 6 = w/2 and 6 = 3n/2. These are the positions of 
maximum and minimum effectiveness and would be 
representative of the range of velocities to be en- 
countered. At6@ = 2/2 the blades were traveling in the 
direction of the wind-tunnel airflow, and at 6 = 37/2 
they were traveling opposite to this flow. Since these 
blades were held rigidly at the hub, they passed 
through a region of high lift and high torque at @ = 
3/2 and of low lift and low torque 180° from this posi- 
tion. 

The average velocity records were taken at 2-in. 
intervals from about r = 4 to 34 in. because there 
seemed to be an interesting reversal of flow past the 
tip of the blades. The plain velocity readings did not 
indicate anything particular, but the direction of flow 
changed so much that attention was drawn to this 
region from the tip outward. 








VELOCITY (FEET PER SECOND) 













12 16 
RADIUS 


Fic. 15. Induced axial velocity vs. radius for horizontal and 
hovering flight at 6 = 1/2. 
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Fic. 16. Velocity under rotor vs. radius for horizontal flight at 
0 = 4/2. 


INDUCED AXIAL VELOCITY (FEET PER SECONn) 


+ 12 24 
RADIUS 





Induced axial velocity vs. radius for horizontal and 


hovering flight at @ = 37/2. 


Fic. 17. 
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Fic. 18. Velocity under rotor vs. radius for horizontal flight at 
0 = 3x/2. 
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From Fig. 15 it can be seen that there is a reversal of 
axial flow from about r = 22 in. outward for the 8.52 
and 13.4 ft. persec.runs. The 19.3 ft. per sec. run does 
not show this negative v’, and this seems to follow the 
trend set by the 13.4 ft. per sec. run, which shows a 
decrease over the slower velocity. Although an axial 
velocity is indicated past the blade tip, not too much 
emphasis should be placed on it, since it was difficult to 
obtain accurate angle readings in this region of dis- 
turbed airflow. These measurements were also taken 
close to the edge of the tunnel, and, although qualita- 
tively they are probably valid, the actual values ob- 
tained are less accurate than the rest. 


Theoretically, the axial flow immediately under the 
rotor for the 19.3 ft. per sec. run should be the lowest, 
and then it should increase as the horizontal velocity 
decreases. The 8.52 ft. per sec. run follows this trend 
as compared to the other two, but the 19.3 ft. per sec. 
run shows about the same axial induced velocity as the 
13.4 ft. per sec. run. It seems that the horizontal air- 
flow is great enough to disrupt this tendency even at 
2 in. below the blades. (This is evidenced again at 
9 = 3/2, as will be explained later.) 

Fig. 16, from which the axial induced velocities were 
taken, shows the maximum effectiveness of the blade to 
be at about 18 in. out from the center or at about 75 
per cent of the radius. There is a gradual rise in in- 
duced velocity up to this peak and then a decline. The 
percentage of increase becomes less as the horizontal 
velocity increases. 

At 0 = 37/2 the effectiveness of the blade should be 
just opposite to that just discussed. The forces on the 
blade here are increased by the horizontal airflow, and 
this increase is supposed to be great enough to not only 
offset the opposite side but to cause the overall effective- 
ness of the rotor to increase. This is evidenced by both 
the force tests and the theory. At all three wind-tun- 
nel velocities a reversal of flow is caused at the tip, the 
values of the 13.8 ft. per sec. velocities shifting between 
the other two. Although theoretically the induced 
axial velocities should increase with the horizontal 
velocity, Fig. 17 shows a decrease. The worst curve 
(19.6 ft. per sec.) could be explained from the fact 
that the blade is evidently stalled. This is indicated 
by Fig. 18, which shows the gradual increase of effec- 
tiveness along the radius for the 13.8 and 8.55 ft. per 
sec. runs but shows no such increase for the other. 
The 19.6 ft. per sec. run indicated a near stalled condi- 
tion, because at times it gave velocity readings that 
varied considerably. The low readings plotted were 
the most consistent. 

“These velocities discussed so far have been average 
values taken under the blades at the positions indicated. 
They ranged from the lowest expected values to the 
highest. Since the force on the blade is a function of 
v’ and since v’ has been shown to be greatest at @ = 
3/2 for any horizontal velocity tried, this fact checks 
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with the theory that predicts this effect. However, 
this theory also only considers an average velocity, and 
it is known that the velocity is not steady, but pulsat- 
ing, because the blades do not form a solid disc. There- 
fore, instantaneous velocity tests were made to try to 
determine the magnitude of these variations. 


Instantaneous Velocities 


Three radial positions at r = 12, 17, and 21 in. were 
chosen, and instantaneous velocity records were made 
at 6 = w/2 and 37/2. Several check points were also 
made at 6 = 0. Several sample records of these runs 
are shown in Fig. 19, the thin horizontal reference line 
through each photograph representing a 15 ft. per sec. 
velocity. These records show the entire Velocity under 
the rotor blades, not just the axial induced velocities. 
Such records were made because no means of deter- 
mining the instantaneous angle was provided. The 
timing wave is 60 cycles per sec., and the blade impulses 
are 22.5 per sec. Therefore, each blade takes up ap- 
proximately three timing wave cycles. 

From these records the impulses were not only dis- 
covered to be large, but it was found that the low veloc- 
ity values of the waves actually went below the hori- 
zontal velocity value in almost all instances. At @ = 
x/2 and with the horizontal velocity at 13.4 ft. per sec., 
the percentage of velocity change based on the hori- 
zontal velocity was at r = 21, 17, and 12 in., respec- 
tively; 23.1, 21.6, and 30.2 per cent. At Y = 19.3, 
the percentages are 13.2, 19.3, and 22 per cent, respec- 
tively. At @ = 32/2 and with the horizontal velocity 
at 13.8 ft. per sec., the percentage change for each simi- 
lar station is 43.1, 34.4, and 36.9 per cent. For 
Y = 19.6 ft. per sec., these values are 19.1, 27.5, and 
33.7 per cent, respectively. These values show that 
the impulses are greatest at @ = 32/2 where the total 
force is also greatest. They also show that the in- 
crease of horizontal velocity causes a definite decrease 
in the percentage change, even though some of the ac- 
tual amounts of velocity change are greater for the 
higher horizontal velocity. 

A close examination of these records showed that at 
6 = 3x/2, where the impulses are the greatest, the im- 
pulse builds up gradually and drops off sharply at 
r = 12 in.; it builds up and drops off at about the 
17 in.; and it has an extremely er- 
ratic behavior at r = 21 in. This is true for both the 
13.8 and 19.6 ft. per sec. runs. At 6 = 2/2 no definite 
impulses like those exhibited on the opposite sides were 
present. In all cases the average velocity measured 
previously with the direction-velocity ‘‘hot wire’’ was 
between the maximum and minimum values obtained 
by the instantaneous wire. 

From this investigation it has been shown that al- 
though the basic theory could be applied to plain hover- 


same rate atr = 
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Fic. 19a (top). @ = 2/2; v = 19.3 ft. per sec.; r = 12 in. 
Fig. 19b (center). 6 = 0; v = 13.8 ft. per sec.; r = 8in. Fig. 
19c (bottom). 0 = 3/2; v = 13.8 ft. per sec.; r = 12 in. 


ing flight, where the design method introduced here 
proves fairly satisfactory, there is much about the 
flow under the rotor during horizontal flight which this’ 
theory does not include. Such effects as the apparent 
loss of v’ at high horizontal velocities would reduce the 
possibilities of obtaining stability by down-wash from 
the rotors. The stalling of the rotor blades during part 
of the cycles and the reversed flow beyond the tips dur- 
ing horizontal flight are two more factors not shown in 
the theory. However, perhaps the most startling find 
is in the instantaneous velocity changes. Although it 
was well understood that impulses were to be expected, 
it is doubtful that the percentages obtained by these 
tests were anticipated. Not only were they much 
greater than seemed evident at first, but the fact that 
they actually went into the negative region, with respect 
to the horizontal flow, shows the possibility of greater 
vibration effects and the need for extra stress considera- 
tion. 
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The Column Paradox 


F. R. SHANLEY* 
Lockheed Aircraft Corporation 


Since this article was received, the author has informed us that he has worked. out the equilibrium equation 


for a simplified, two-flange column. 
with increasing lateral deflection. 
to the value predicted by the reduced modulus theory. 


preparation.—EDITOoR. 


INTRODUCTION 


AS SHOWN BY VON KArRMAN in his classical paper,' 
column action in the inelastic range (beyond the 
proportional limit) involves two different moduli of 
elasticity: the elastic (Young’s) modulus, and the 
tangent modulus (local slope of stress-strain diagram). 
The “reduced modulus”’ is a function of both of these 
and is also affected by the shape of the cross section. 
By using the reduced modulus in place of Young's 
modulus, Euler’s equation has been generalized to in- 
clude stresses well beyond the elastic limit. 


REDUCED MopuULUS THEORY 


If the tangent modulus is used directly in the Euler 
formula, the resulting critical load is somewhat lower 
than that given by the reduced modulus theory. This 
simpler formula, originally proposed by Engesser, is 
now widely used by engineers, since it gives values that 
agree very well with test data. Nevertheless, the re- 
duced modulus formula is still generally considered to 
be the true theoretical solution for perfect columns, 
and the lowér test values are explained by references 
to unavoidable eccentricities, testing technique, and 
other errors. 

The reduced modulus theory is derived by assuming 
that, after the column reaches the critical uniform stress, 
the column bends. This causes the strain to decrease 
on one side and increase on the other. For the in- 
creasing strains the resulting stresses are given by the 
tangent modulus, while for the decreasing strains the 
elastic modulus gives the relation between strain and 
stress. Tests have shown that, when the strain reverses, 
actual conditions are in accordance with this theory. 


DISCUSSION 


But there is an implied assumption in the derivation 
of the reduced modulus theory that is open to queStion. 
It is, in effect, assumed that something keeps the 
column straight while the strain increases from that 
predicted by the tangent modulus theory to the higher 
value derived from the reduced modulus theory. 
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The column starts to bend at the tangent modulus load and the load increases 
If the tangent modulus is assumed to be constant, the column load is asymptotic 
This development will be covered by a second paper, now in 


Actually, there is nothing (except the column's bending 
stiffness) to prevent the column from bending simul- 
taneously with increasing axial loading. Under such 
conditions the compressive strain could increase on one 
edge of the column while remaining constant on the 
other, or it could increase at a different rate on opposite 
edges. If such action were assumed, the tangent 
modulus would apply over the entire cross section, and 
the theoretical buckling load would be that predicted 
by the tangent modulus theory. This creates a paradox, 
because, if all of the strains equal or exceed the tangent 
modulus value, the average strain will be greater than 
that predicted by the tangent modulus theory. 

The assumptions involved in the reduced modulus 
theory also represent a paradox. The theory predicts 
that the column will remain straight up to the calcu- 
lated maximum load, but it also shows that some 
strain reversal is needed in order to provide the addi- 
tional column stiffness required beyond the tangent 
modulus load. It is impossible to have strain reversal 
in a straight column. 

Evidently, the maximum column load will be 
reached somewhere between the loads predicted by the 
two theories. The entire problem should be reviewed 
on the basis that axial loading and bending can occur 
simultaneously. The use of the principle of super- 
position is not valid in this case. 


CONCLUSIONS 


Since it will require some reversal of strain in order 
to exceed the tangent modulus load, it can be predicted 
that bending will begin as soon as that load is exceeded. 
The tangent modulus theory can therefore be regarded 
as correct for predicting the maximum axial load at 
which a perfect centrally loaded column will remain 
straight. 

To find the maximum load for a perfect column in 
the plastic range, it will evidently be necessary to 
analyze each case separately, using a step-by-step 
process and employing the general methods followed 
by von Karman in his solution of the eccentric column. 
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Errata 


Tailless Airplanes and Their Practical Application, 
THE AERONAUTICAL SCIENCES on pages 525-536. 


(1) Page 530, Eq. (8 


(3) Page 534, Eq. (37) 
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(4) Page 535, under the heading ‘‘Conclusions,”’ 





We regret the following errors in the article by Maurice A. Garbell, 
” which appeared in the October, 1946, issue of the JouURNAL OF 


dC 
), the left-hand expression should be: ( *). 


. or sweptback high-aspect-ratio wings and in tail-first airplanes; 


“Theoretical Principles of Wing-Tip Fins for 


at 


° IC 
(2) Page 531, Eq. (12), the left-hand expression should be: — 
d(rb/2 2V) ) 


, the first member of the bracketed expression, right-hand side, should be: 
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the 8th and 9th lines should read: 


it is maintained that the central vertical 














